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Summary—Fine particles of graphite (90 to 720 «) were burnt at 500-600°C. At the same 
; time the rate of reaction and the B.E.T. area of the sample were determined. While the experi- 
953 ment is taking place the rate of combustion per cm? of B.E.T. surface keeps rising. The evolution 
of the fine porosity of the particle allows an interpretation of this variation. It is only at the 
end of the combustion when the particle has become very cavernous that the measured rate 
of reaction has a definite physical significance. The complications which can happen according 
to the nature of the porosity and its distribution between the particles are examined. This 
method and this interpretation are applicable to any gas-solid reaction of the carbon combustion 
type; they particularly deal with the problem of coke reactivity. 


Résumé—Les auteurs décrivent leur méthode qui consiste 4 braler 4 basse température (500- 
600°C) jusqu’a disparition compléte un échantillon de graphite en trés petits grains (90 a 720 yu) 
calibrés, en déterminant parallélement la vitesse de réaction et la surface B.E.T. de l’échantillon. 
Ils donnent pour exemple une courbe de combustion qui correspond 4 un cas relativement simple. 
La vitesse de combustion rapportée au cm? de surface B.E.T. ne cesse d’augmenter au cours de 
l'expérience. Ils montrent comment l’évolution de la porosité fine du grain permet d’interpréter 
cette variation. C'est seulement en fin de combustion, lorsque le grain est devenu trés caverneux 
qu'on peut atteindre et mesurer une vitesse de réaction qui a une signification physiquement 
définie. Ils examinent rapidement les diverses causes de complications qui peuvent survenir 
suivant la nature de la porosité et sa distribution entre l'ensemble des grains. La méthode et 
l'interprétation intéressent toute réaction gaz-solide du méme type que celle de la combustion du 
carbone ; en particulier, elles intéressent le probléme de la réactivité des cokes. 


Les valeurs données pour la vitesse de combustion 
du carbone, la composition des produits formés 
ainsi d’ailleurs que les caractéristiques cinétiques 
de la réaction varient beaucoup d'un échantillon 
a autre. C’est ainsi que J. R. Artuur [1] décrit 
la combustion de deux lots de carbone (un graph- 
ite artificiel et un coke) qui brdilent a la méme 
vitesse alors que l'un est vers 900°C et l'autre vers 
600°C ce qui, pour une énergie d’activation 
d’ArrueEntvus de 40 K. cal. correspondrait, toutes 
choses égales d’ailleurs, A un rapport de 1 a plus 


9 


de 10° entre les constantes de vitesse attachées, 
pour la méme température, & ces échantillons 
respectivement. 

Nous avons donné [2] de premiéres raisons de 
penser qu’indépendamment de la structure cris- 
tallographique (graphite ou carbones amorphes) 
et de la pureté du carbone, la porosité de la masse 
soumise & la combustion joue un réle important 
dans cette dispersion des résultats. La présente 
étude confirme cette opinion. Elle décrit notre 
méthode d’étude et propose une interprétation de 
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influence de la porosité fine du combustible sur 
la vitesse de combustion et sur l’évolution de 
cette vitesse au cours de l’expérience. 

La méthode expérimentale a été deécrite 
antérieurement [3]. Elle se 
conditions habituelles d'utilisation des combusti- 
bles solides a la différence qu’elle est & trés petite 
échelle afin que les conditions expérimentales, 
température, soient strictement 


rapproche des 


pression et 
définies. 


L’appareil est en somme un trés petit fourneau 
construit en verre pyrex ou en silice fondue. Un 
courant d’air sec et décarbonaté traverse a la 
pression ordinaire et A un débit de 5 A 15 litres a 
lheure une pastille de silice frittée qui porte une 
couche (de 3 A 6 mm de hauteur et de 2 A 3 cm 
de diamétre) de grains de graphite Acheson de 
haute pureté. 

Ces grains sont préparés par broyage au mortier 
d’agate de baguettes, barres plaques de 
graphite,* et séparés par tamisage. 
généralement employé deux calibres de grains le 
diamétre moyen étant 720 u (passant au tamis 
A.S.T. M.20 et refusé au tamis 30) et 90 » (passant 
au tamis 140 et retenu sur tamis 200). Les 
volumes de ces grains sont done dans le rapport 
de 500 & 1 environ. 

La température de la charge (fixée A + 0-5° 
prés) est choisie pour que la combustion soit trés 
lente (plusieurs dizaines d’heures pour briler 
0-5 gr.). Cette condition est indispensable, comme 
on l’a reconnu depuis longtemps, pour étre assuré 
que la vitesse mesurée n’est pas celle de phéno- 
ménes occasionnels de diffusion du réactif gazeux 
dans le lit de combustible. Nous allons montrer 
toutefois, ce qui semble-t-il est original, que cette 
condition n’exclut pas le réle, qui est déterminant, 
de la diffusion a l’intérieur du grain. 

Pour ces basses températures, le gaz qui sort 
du four contient un large volume d’oxygéne qui 
n’a pas réagi, environ 1% en volume de CO? et 
0-1 & 0-4% en volume de CO. C’est dire qu’il se 
fait & la fois CO? et CO. Le rapport des con- 
centrations de ces oxydes varie beaucoup d’un 


ou 
Nous avons 


* Nous adressons nos remerciements & la Compagnie 
PECHINEY qui a bien voulu nous fournir gracieusement ces 
échantillons., 


graphite a l'autre; par exemple, le rapport 
[CO] 

[CO] + [CO,] 

que nous avons longuement etudié alors qu’il est 


est 0-4 pour un graphite artificiel 


inférieur & 0-1 pour un graphite naturel (Madas- 
gascar) dans les mémes conditions. Nous n’avons 
pas encore examiné la question de savoir si cette 
composition correspond aux produits primaires 
de la réaction, ce qui parait d’ailleurs peu pro- 
bable. 

La teneur en CO et CO? du gaz de réaction (que 
nous déterminions antérieurement par micro- 
gravimétrie [3] est mesurée maintenant de facgon 
continue par absorption du rayonnement i.r. Il 
est facile d’en déduire la carbone 
gazéifiée en un temps donné (une min parexemple); 
cette masse définit la vitesse de combustion V. 


masse de 


La figure 1 donne cette vitesse en fonction de 
la durée de lexpérience ; elle représente A titre 
d’exemple la combustion a 609°C de 0-5 gr. d'un 
échantillon (immatriculé AFC.1) graphite 
Acheson pur (0-02% de cendres) sous les deux 
calibres de grains. Ces grains n’avaient subi aucun 
traitement (de dégazage par exemple) avant 
d’étre soumis & la combustion. II est important 
de remarquer sur les courbes de cette figure que 
lexpérience a été poursuivie jusqu’é disparition 
compléte du combustible (V = 0); nous verrons 
bientét l’intérét de cette fagon de procéder. 


de 
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D’autres échantillons, également de haute 
pureté et d’origine analogue, ont été essayés dans 
les mémes conditions. Les courbes obtenues bien 
qu’elles présentent des analogies sont nettement 
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différentes entre elles, en particulier du point de 
vue quantitatif. En fait, chaque lot eramind fournit 
une courbe qui lui est particuliére. 

D’autre part, pour un lot déterminé, on est bien 
embarrassé pour définer une vitesse de combustion 
caractéristique puisque cette vitesse ne cesse de 
varier dans d’importantes proportions d'un bout 
& l'autre de l’expérience comme le montre la 
figure 1. Ainsi, les combustions complétes aux- 
quelles nous nous sommes astreints montrent 
linanité des vitesses stationnaires’’ (en fonc- 
tion du temps) qu’on pourrait espérer atteindre 
par expérience directe. 

De fait, pour une réaction hétérogéne, |'expres- 
sion d’une vitesse de réaction n’a de sens que si 
elle est rapportée a l’aire qui participe a la 
réaction. Afin de connaitre celle-ci, la surface 
réelle de la charge de graphite a été mesurée grace 
& la méthode classique de BrunaverR, Emmett 
et Te_Ler (B.E.T.) par adsorption de méthane 
(ou d’azote) & — 195°C.* 

La mesure effectuée sur la charge de graphite 
avant combustion montre que l’aire réelle des 
grains est beaucoup plus grande que leur surface 
géométrique; pour le lot AFC 1 auquel se 
rapporte la fig. 1, le “ facteur de rugosité”’ est 
70 sur le grain fin et 220 environ sur le gros 
grain. 

Le grain est done poreux. Cette observation 
essentielle est recoupée par d'autres faits : 
branches d’hystérése que manifestent les iso- 
thermes d’adsorption dans la mesure des surfaces 
B.E.T., lenteur de pénétration du benzéne dans 
le grain lors de la détermination picnométrique 
de densité du graphite, diffusion aux petits angles 
des rayons X par le graphite. 

Aprés qu'elle a subi une combustion d'une 
durée déterminée, la charge de graphite est sortie 
du four, pesée et on en mesure A nouveau sa 
surface B.E.T. De la sorte, la construction d’une 
courbe compléte de combustion exige une longue 


* Ces mesures ont été effectuées par Mr. Lucien 
BoNNETAIN que nous remercions bien vivement. Signalons 
& ce propos que le type nouveau et trés particulier d’iso- 
thermes d'adsorption qui a été découvert a cette occasion 
(L. Bonnetain, X. Duvat et M. Lerort C.R. 1952 234, 
1363) vient d’étre confirmé indépendamment (M. H. 
Pottey, W. D. Scuarrrer et W. R. Smitru, J. Phys. 
Chem., 1953 57 469). 


expérimentation, mais il parait indispensable de 
consentir cet effort si l'on veut comprendre le 
mécanisme de la combustion et mesurer une 
vitesse qui soit physiquement définie. 


Le document le plus important que fournit 
cette étude est représenté sur la figure 2 ot la 
vitesse de combustion V’ (en mg de carbone 
gazéifié par min) du lot AFC 1 est portée en 
ordonnée et la surface B.E.T. (m?*) en abcisse S. 
Pour étre complet ce graphique pourrait com- 
porter suivant une troisiéme coordonnée I’indica- 
tion de la masse brilée & chaque instant : celle-ci 
est représentée par les diverses valeurs m portées 
sur le plan et qui varient de m = 0-5 gr (instant 
initial) A m = 0. 


On distingue aisément sur ces courbes trois 
parties que nous commenterons pour exemple 
sur la courbe en trait plein qui représente la 
combustion du grain 20-30. Cette courbe 
montre : 

1°/ Une phase initiale qui s’étend jusqu’A ce 
que 20°, de la charge soit brilée (m = 0-1 gr); la 
surface B.E.T. ainsi que la vitesse de réaction 
augmentent et, approximativement, dans le méme 
rapport. 

2°/ Une phase intermédiaire qui va jusqu’- 
environ 80%, de la masse brilée (m = 0-4 gr); S$ 
diminue considérablement tandis que V varie 
trés peu. II faut spécialement remarquer au début 
de cette phase une partie de la courbe qui apparait 
peu sur cette figure bien qu’elle corresponde a 
20% de la diminution de masse (de m = 0-1 gra 
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m = 0-2 gr) et pour laquelle S commence a dimi- 
nuer alors que V continue &a croitre. 


3°/ Une phase finale qui se termine nécessaire- 
ment a l’origine des coordonnées; S et V dimi- 
nuent concomitamment, la porportionnalité de 
ces deux grandeurs étant exactement vérifiée dans 
la partie extréme de la courbe, c’est-A-dire en fin 
d’expérience. 


Il est facile de reprendre les points expéri- 
mentaux de cette figure sous une autre forme. 
Dans ce nouveau graphique, (fig. 3) nous donnons, 
en fonction de la masse consumée, la vitesse de 


réaction V° rapportée au m* de surface B.E.T. 


mg/min/m? | 609°C) 


20.30 


140.200 


VITESSE 
COMBUSTION 


° 


° o2 o4gr Os 
MASSE BRULEE 


Fig. 3. 


Nous retrouvons évidemment les trois phases que 
nous venons de dire. Jusqu’environ 20% de 
masse brilée, cette vitesse V° est & peu prés 
constante et de l’ordre de 0-2 mg/min. m*; puis 
elle augmente assez vite pour se fixer a la valeur 
constante de 1-5 mg/min qui n'est d’ailleurs 
atteinte que lorsque plus de 95% de la charge 
sont disparus. I] est intéressant d’observer que 
cette valeur limite est la méme pour les deux 
calibres de grains. 


L’essentiel de |'interprétation tient dans l’idée 
trés simple que toute l'aire mesurée par la méthode 
B.E.T. ne participe pas nécessairement a la com- 
bustion. 

Sans doute, toute cette surface est-elle accessible 
aux gaz puisque la mesure par adsorption qu’on 
en fait est précisément fondée sur l’accés d'un 
gaz a la surface, et il n'y a aucune raison de penser 
que la surface décelée par adsorption ne soit pas 
également accessible aux gaz & la température 


de combustion. Encore faut-il pour que la com- 
bustion ait lieu en un endroit donné de la surface 
que l’oxygéne puisse parvenir jusqu’aé cet endroit. 
Or, l’oxygéne peut étre complétement consommé 
avant de l’atteindre. 

D’une facon plus générale, suivant |'endroit 
considéré de la surface, les conditions expéri- 
mentales & un instant donné de la combustion 
peuvent étre trés différentes. A ce propos, nous 
avons A examiner trois facteurs: nature (et 
état) de la surface a l’endroit considéré, tempéra- 
ture et pression de l’oxygéne en ce lieu. 

On ne voit pas de raison pour que l’activité 
propre de la surface ne soit pas la méme en chaque 
point; nous admettrons donc que la fraction 
active (qui peut étre trés petite) de l'unité d’aire 
est partout identique a un instant donné (ce qui 
ne veut pas dire, que celle-ci reste nécessairement 
égale & elle-méme pendant toute la durée de 
l’expérience). 

Nous négligerons de méme l'éventuelle différ- 
ence de température qui pourrait exister entre le 
coeur et la périphérie du grain du fait de l’exo- 
thermicité de la réaction, le centre pouvant étre 
plus chaud que la surface. Nous avons toutes 
raisons de croire en effet qu’A ces trés faibles 
vitesses de combustion, l’équilibre thermique est 
réalisé entre tous les grains de la charge et dans 
chaque grain considéré isolément. 

Par contre, si lon veut bien considérer que 
loxygéne qui diffuse depuis l’entrée d’un pore 
vers le fond (tandis que les produits de la réaction 
diffusent en sens inverse) est nécessairement 
consommeé pendant ce cheminement, on doit 
admettre que, suivant la profondeur du point 
considéré sur la paroi du pore, la pression de 
l'oxygéne peut prendre toutes les valeurs depuis 
la pression de ce gaz dans l’air a la pression 
ordinaire (1/5 atm) jusque zéro. 

La vitesse globale qu’on mesure étant la somme 
de toutes les vitesses locales, il conviendrait de 
connaitre en chaque point de la surface la valeur 
de la vitesse de combustion en fonction de la 
pression d’oxygéne en ce point. A ce titre, nous 
nous représentons la combustion du grain de la 
facon suivante. 

Jusqu’é une certaine profondeur dans le pore, 
la pression d’oxygéne reste suffisante pour que 
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l'acte d’adsorption de ce gaz soit extrémement 
rapide comparativement a la vitesse de trans- 
formation chimique en phase adsorbée ou a la 
vitesse de désorption. Jusqu’a cette valeur de la 
pression c'est & dire jusqu’a cette profondeur dans 
le pore, la vitesse de réaction serait done indeé- 
pendante de la pression d’oxygéne (ordre zéro) et 
identique en tout point de la paroi. II existerait 
ainsi autour du grain une sorte décorce, plus ou 
moins épaisse suivant la température, & l'intérieur 
de laquelle toute la surface serait soumise & une 
vitesse bien définie de combustion. Cette écorce 
est vraisemblablement d’autant plus mince que 
la température est plus élevée puisque la con- 
sommation de l’oxygéne qui diffuse dans un pore 
est évidemment d’autant plus rapide que la 
température est plus grande. 

Au dela de cette écorce en allant vers l'intérieur 
du grain, la pression d’oxygéne se trouvant trés 
sensiblement diminuée, l’ordre zéro ne serait plus 
respecté mais tendrait progressivement, & mesure 
que la pression ps2 locale s’abaisse, vers lordre 
1 qui serait effectivement atteint pour les plus 
basses pressions d’oxygéne. Ainsi, sous l’écorce a 
Vintérieur de laquelle chaque cm? de surface 
briile A la méme vitesse, il existerait une zone ot 
la vitesse locale de combustion diminuerait 
progressivement jusqu’é devenir nulle.* 

Ainsi la vitesse globale mesurée (qui représente, 
nous le rappelons, une masse de carbone gazéifié 
par min) serait la somme de deux termes que nous 
désignerons pour simplifier par vitesse de com- 


‘ 


‘en régime externe ”’ et vitesse de com- 


“ 


bustion 
bustion “en régime interne,’ ce second terme 
¢tant lintégrale définie d'une infinité de vitesses 
locales variables suivant le point considéré. La 
vitesse locale moyenne de la zone interne étant 
nécessairement plus petite que la vitesse locale 
(identique en chaque point) de la zone externe, 
pour que ce second terme ne soit pas négligeable 
vis & vis du premier, il faut penser que la zone 
de combustion, “ en régime interne ”’ est épaisse 


* Le calcul permet de préciser qu’é mesure qu’on se 
déplace depuis Pentrée du pore vers le fond, poz décroit 
d’abord suivant une parabole dans toute la zone ou l'ordre 
zéro est vérifié ; puis cette courbe tend vers la forme d'une 
exponentielle qui serait effectivement atteinte lorsque la 
réuction devient d’ordre 1. 
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par rapport & l’écorce qui brile 


en régime 
externe.”’ 

Cette représentation permet de rendre compte 
trés simplement des résultats expérimentaux. 

Ainsi, pour la courbe que nous avons commentée 
(fig. 2: grain 20-30), dans la phase initiale de la 
combustion, il s’exerce une intense corrosion du 
grain qui multiple par 3 environ la surface de 
celui-ci (pour 0-5 gr de graphite, S passe de 0-5 
i 1-4m?). La réaction se fait essentiellement en 
La propor- 
tionnalité entre V et S est approximativement 


régime de “ combustion externe.” 


respectée, surtout en début de combustion. 

Mais le développement de la surface ne se 
poursuit pas indéfiniment. L’enlévement de 
matiére intervient en sens inverse de l’effet de 
corrosion : le détail du phénoméne est évidem- 
ment fort complexe mais, au total, on constate 
qu’entre 20 et 80% de la consommation com- 
pléte de la charge, la surface diminue de 1-4 a 
0.35 m*, Cependant la vitesse de combustion ne 
varie que trés peu pour la raison que la corrosion 
a augmenté et ne cesse d’augmenter I’accessibilité 
des pores & loxygéne. De facon générale, les 
parties trés enfoncées de la surface qui ne brilaient 
pas, tendent a mesure qu’elles se découvrent, a se 
consumer d’abord en régime de combustion 
interne,”’ puis, le pore s’étant creusé davantage, 
‘externe c’est a dire 


‘ 


en régime de combustion 
au maximum de vitesse. C'est la le fait essentiel. 

Enfin, le grain étant devenu trés caverneux, 
toute la surface brile au maximum de vitesse. La 
pression locale poz est la méme en chaque point. II 
y aurait-il des différences de pression d’un endroit 
& un autre que celles-ci ne se traduiraient plus sur 
la vitesse locale qui obéit partout a l’ordre zéro. 
On est ici dans la phase finale de la combustion. 
Elle correspond & un cas trés simple. En effet, 
tout se passe comme si on substituait au grain 
caverneux un grain absolument compact et de 
méme surface (réelle). La proportionnalité entre 
V et S est exactement vérifi¢e pour les derniers 
temps de I’expérience. 

Au total, la corrosion du grain permet done 
d’interpréter de fagon rationnelle |évolution de 
la vitesse de combustion en fonction du temps. 
Dans |’état actuel de nos recherches, rien n’auto- 
rise cependant A croire que la corrosion du grain 
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soit la raison exclusive de cette évolution. On 
peut concevoir que, parallélement, l'état de 
surface de l’unité d’aire se modifie spontanément, 
par exemple par formation d'une couche adsorbée 
inhibitrice. Des recherches dans ce sens ne 
pouvaient toutefois étre .entreprises utilement 
avant que l'influence de la porosité n’ait été 
elle méme clairement reconnue. 

Il est important de remarquer enfin que c’est 
seulement dans la derniére partie de la phase 
finale de la combustion que la vitesse observée 
prend une signification physiquement deéfinie car, 
& ce moment, la surface B.E.T. participe totale- 
ment et de la méme facon & la combustion. Pour 
les autres parties de la courbe de combustion, la 
distribution des vitesses locales de réaction et la 
part de la surface B.E.T. qui intervient réellement 
dans la combustion, sont inconnues. 


On arrive ainsi & la conclusion d’apparence 
paradoxale (du moins pour le graphite que nous 
étudions ici ou des graphites analogues) que c’est 
seulement lorsque la quasi-totalité de la charge 
est brailée qu'il est possible d’atteindre et de 
mesurer la vitesse de l'acte chimique dans des 
conditions définies, grandeur que nous proposons 


d’appeler “ vitesse réelle de combustion” V,, 


Cette conclusion justifie & posteriori notre méthode 
qui consiste A briler toute la charge. 


Pour le graphite AFC.1 que nous venons 
d’examiner, cette vitesse absolue est égale a 
1-25 10'* atomes,cm? sec A 609°C. 

Toutefois, il semble qu’il subsiste des différences 
sensibles entre les vitesses réelles de combustion 
V, de divers graphites. L'examen systématique 
de nombreux loig de provenances diverses per- 
mettra de préciser ce point. Ces différences ne 
doivent pas surprendre ; il est concevable en effet 
que la fraction active de l'unité de surface dépende 
de la structure cristallographique du carbone 
(degré de graphitation). L’état de surface en fin 
d’expérience peut d'autre part varier d'un 
graphite a l'autre. 

Quoi qu'il en soit, ces écarts (qui, dans |'état 
actuel de nos recherches, paraissent étre de l'ordre 
d'un facteur de 2 ou 3) sont petits comparés aux 
différences beaucoup plus considérables que l'on 
constate entre des vitesses brutes de combustion 


et que nous attribuons essentiellement aux effets 
de la porosité. 

L’exemple du graphite AFC. A a été choisi car 
il correspond a un cas relativement simple et 
fondamental. D’autres échantillons fournissent 
des résultats plus compliqués et dont l’interpré- 
tation est moins simple. Toutefois, l’analyse qui 
précéde du mécanisme de combustion du grain 
reste non seulement valable mais elle constitue 
indispensable clef des nouveaux problémes qui 
se posent. 


En effet, comme on peut logiquement l’attendre 
au terme de cette analyse, les différences de 
porosité doivent nécessairement se traduire dans 
les courbes de combustion. Ces différences — et 
les complications qu’elles entrainent — paraissent 
provenir essentiellement de trois causes : 


1°/ Le diamétre initial moyen des pores est 
plus ou moins large, l’influence de la diffusion dans 
les pores étant d’autant plus marquée que la 
porosité est plus fine. 

2°/ Les pores ont grossi¢rement le méme dia- 
métre ou, au contraire, les dimensions de ces 
diamétres sont largement dispersées. 

3°/ Les pores sont distribués de facon sensible- 
ment uniforme dans la masse du graphite ou, 
au contraire, certains grains sont notablement plus 
poreux que d'autres, de sorte que la charge se 
comporte comme un mélange de diverses qualités 
de grains. Cette hétérogénéité de distribution 
entre les grains intéresse le nombre des pores et 
(ou) leurs dimensions.* 

Le probléme peut done se présenter sous des 
aspects extrémement complexes. Nous n’exami- 
nerons pas ici ces complications. C'est essentielle- 
ment la méthode de travail que nous désirions 
exposer. Nous indiquerons seulement que le 
graphite (AFC.1) dont nous avons commenté la 
combustion posséde des pores trés fins, & peu prés 
analogues entre eux et distribués de facon uni- 
forme dans la masse; ces deux derniéres carac- 
téristiques justifient que nous le considérions 
comme un cas relativement simple. 


* De fuil, des diverses causes de complications ne sont 
pas indépendantes ; ainsi, il semble que la présence de 
pores trés larges s accompagne toujours d'une distribution 
non uniforme de la porosité entre l'ensemble des grains. 
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Il va sans dire enfin que cette méthode de 
travail et cette interprétation ne se limitent pas 
au cas graphite que nous n’avons choisi que pour 
disposer d’un combustible bien défini et éliminer 
ainsi (ou du moins considérablement diminuer) 
linfluence d'autres facteurs tels que la structure 
cristallographique et la pureté chimique. Elles 
intéressent toutes les réactions hétérogénes du 
méme type (reaction gaz — solide avec formation 


de produits gazeux). En particulier, dans le 
domaine de la combustion du carbone, elles 
intéressent directement une étude rationnelle du 
probléme depuis longtemps posé de la réactivité 
des cokes. 

L’un de nous (R.M.) adresse ses trés sincéres 
remerciements a la Cie Pechiney pour l’assistance 
financiére qu’il a recu de celle-ci pendant toute 
lexécution de ce travail. 
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Mixing of solid particles 


R. Buumperc and J. S. Marirz* 
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Summary — The theoretical concepts of complete and incomplete mixing have been discussed, 
and a definition for complete mixing has been proposed, using a statistical hypothesis as the test 
criterion. The application of this has been confirmed experimentally and it has been shown 
that, with random sampling, n, the number of particles in a sample, and z, the proportion of 
particles of type A in the sample, follow normal distributions. 

It has further been shown that a small number of samples examined statistically, is sufficient 
to indicate whether or not mixing is complete at any time during the mixing process, and that 
the progress of mixing with time may be followed statistically. 


Résumé—Aprés avoir discuté les concepts théoriques de mélange imparfait et de mélange 
parfait,” l'auteur propose une définition du “ mélange parfait posant @ la base du critére 
expérimental une hypothése statistique. L'expérience est en accord avec cette hypothése : elle 
montre que le nombre de particules n, est la proportion z des particules du type A dans l’échan- 
tillon, obéissent A une loi de distributions normales, lorsque les échantillons sont choisis au 
hasard. 

De plus, l'auteur établit que, durant l'opération du mélange, l'examen statistique d'un petit 
nombre d’échantillons indique si le mélange est ou n'est pas parfait & un instant donné : il est 
done possible, par l'observation statistique, de suivre les progrés du mélange en fonction du 


temps. 


1. INTRODUCTION 

Although during the past 20 years a large number 
of papers have been published dealing with various 
aspects of mixing, the fundamentals of solid-solid 
mixing have barely been touched. In view of the 
application of the dry mixing operation to many 
widely different processes, it was felt that addi- 
tional fundamental work was warranted. 

A study of the literature showed that there 
have been two lines of approach to the problem 
of solid-solid mixing thus far. Brotruman and 
co-workers [1] and Coutson and Marrra [2], [3] 
approached the problem from the kinetics view- 
point, while Lacey [4], Busiix [5] and De Jone 
[6], all adopted the statistical approach using 
the value of the standard deviation as the mixing 
criterion. Of all these workers, only CouLson and 
Mairra, and Lacey, published experimental 
results in support of their hypotheses. 

Lacey was concerned with degree of mixing, 
and in his experimental work used small batches 


of relatively large (0-2in. diameter) coloured 
spheres. After mixing, he examined the whole 
batch by dividing it into small samples, without 
disturbing the relative positions of the particles, 
and counted the particles in each sample. Since 
the entire batch was examined, the question of 
sampling did not come up for consideration. Lacey 
concluded that it was probably impossible to 
express the state of a mixture adequately by a 
single number, but that it could be expressed by 
a line on a suitable graph, and that this would 
provide a quantitative as well as a qualitative 
method of comparing mixtures and the _per- 
formances of mixing machines. 

CouLson and Mairra examined the main fac- 
tors influencing the rate of mixing and concluded 
that their basic kinetic mixing equation applied 
over a wide range of conditions, but had certain 
definite limitations. 

The kinetic approach offers considerable possi- 
bilities of practical application, but CouLson and 


* rresent address ; Transvaal Chamber of Mines Research Laboratories, Kew Road, Richmond, Johannesburg. 
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Marra omitted to refer to the statisical standards 
on which their conclusions were based; they 
simply state that, if n samples out of 30 were 
found to have approximately the same composi- 
tion as the whole system, then the mixture was 
n 

80 

be applied in practice, because of the impossibility 


< 100%, mixed. Lacey's method cannot 


of sampling and counting an entire batch. 

It was felt that the first stage in a fundamental 
examination of the unit operation of solid-solid 
mixing should be to define complete mixing, to 
settle the question of adequate sampling, and to 
arrive at a means of expressing to what extent 
a mixture approached the completely mixed 
state. 

2. EQuipMENT 

The purpose of these experiments being an 
investigation of mixing fundamentals, and not of 
mixer characteristics, it was considered advisable 
to use a very simple type of mixer. The design 
of Coutson and Marrra was followed in this 
connection, a simple Perspex drum 6 in. diameter 
x 9in. being used. The assembled equipment is 
shown in Fig. 1. 

For the same reason it was thought desirable 
that differences in physical properties of the 
materials being mixed should not enter into the 
investigation at this stage. A closely sized batch 
of washed Cape sand of uniform quality was there- 
fore divided into two portions, the one dyed blue 
and the other red ; these coloured sands were used 
for the experiments. They differed in colour 
only, and this property could have no effect on 
the mixing of the particles. A summary of 
physical properties of the sand is given in Table 1. 


Table 1. Physical properties of Cape sand. 


Particle size range : Tyler mesh 35 +48)—-48 +65 

Bulk density : g/cc 1-60 1-46 

Specific gravity : 20°C | 2-631 2-604 
39 44 


el 
%, Voids 


The sampling device is shown in Fig. 2. It 
consisted of a rod 0-132in. diameter with 9 
cavities of 0-0635 in. diameter drilled at 1 in. 
intervals along its length to a depth of 0-0635 in. 


The first cavity was 0-43 in. from the end. The 
rod was a sliding fit inside an outer sleeve, which 
had a strip 0-076 in. wide cut away along its 


= © 


Fig. 2. Sampler. 


length so that the cavities in the rod could be 
either covered or open. The sampler guide, 
shown in Fig. 8, could be fitted across the top 


Fig. 3. Sampling guide. 


of the drum along any diameter, and the sampler 
inserted. The fact that the material being mixed 
consisted of discrete particles created distinct 
sampling problems and necessitated a standard- 
ised sampling procedure throughout, so as to 
cause as little disturbance as possible in the 
mixture. 
3. CompLeTe MIXING 


Fundamental ideas concerning the completely mixed 
state 

Before it was possible to examine the progress 

of mixing or the degree of mixing achieved, it 

was considered necessary to define the com- 

pletely mixed state. 
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For the experimental work, the sampling 
instrument was so designed that it withdrew at 
each position a very small proportion of the total 
batch, and the sampling could therefore be 
regarded as having taken place at a “point” in 
the batch. The first assumption which was made 
was that if, after a certain set of operations, a 
sample were withdrawn consisting of n particles 
of which r were red, and if the same procedure 
were repeated a number of times, in each case 
starting from the same datum level, then for a 
fixed value of n, corresponding to each sampling 
point, r would be a statistical variable with 
distribution 


p(r) (1) 


@ being a parameter called “ the probability of 
drawing red at a certain point.” 

In the experimental practice it was impossible 
to draw exactly » particles at each point, so that 
n also was a statistical variable. The joint dis- 
tribution of r and n could therefore be expressed 
by 

n 
p(r,n) = (joa p(n) (2 

Making use of the normal approximation to 

the binomial distribution [7], and writing 
a= —= the proportion of red particles inthe sample. 


n 
the joint distribution of z and n became 


nd 
V2 — 8) 
—@)) 


exp 20(1 8) p(n)drdn (3) 


p (a, n) dadn = 


In the derivation of this distribution it was 
assumed that the distributions of both r and n 
could be approximately represented by continuous 
ones. 

From equation (3) the following results were 
obtained : 


E (x) = p, (2) = @ (as was to be expected) 
Hy = — @) E(n") 
Hg (2) 0 
(7) = 367 (1 — E(n?) 


consequently 
B, (2) =0 
36? (1 0)? E(n*) 
B. 


3E (uw?) 
[E 
By making the reasonable assumption that 


var (nn) 
would be small, it could be shown that 


[ 
1 
1 
E(n? 
from which 
B, (2) 3 


i.e., the distribution of 2 would be approximately 
normal with 
mean 0 


—@) 
E (n) 


Variance 


It seemed reasonable to define the completely 
mixed state as that which would obtain when the 
probability of drawing red at a certain point was 
the same at all points in the mixer, Le., for a 
completely mixed batch, the values a7, ....2, 
obtained by random sampling from the batch, 
could be regarded as k independent statistical 
variables each with the same normal distribution. 

The first set of experiments was designed to 
obtain some empirical verification of the above 
theories. 


EXPERIMENTAL PROCEDURE 


The drum was held in the vertical position and 
red sand, 500 g, poured in evenly on the bottom. 
Blue sand, 500 g, was then poured in a layer of 
uniform thickness on top, and the drum attached 
to the stationary turntable. The angle of tilt 
of the drum and the speed of rotation were 
arbitrarily selected, taking into account the 
In all this 
work the drum was rotated at 55 r.p.m. at an angle 


findings of Coutson and Marrra. 


of 30° from the horizontal until it was judged 
that the mixing was complete. The mixer was 
then stopped, and random samples were extracted 
using the sampler. The individual samples were 
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Mixer and auxiliary equipment 
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transferred to separate cavity slides, the particles 
projected on to the s¢reen of a Vickers projection 
microscope, using oblique incident light and 
magnification 12 x, and the number of red and 
blue particles in each sample counted. The small 
quantity of sampled sand was returned to the 
drum, rotation started and continued for a 
similar period of time. Then the mixture was 
again sampled. In this type of test at least 160 
samples were examined and the particles counted. 


Discussion oF RESULTS 
When a distribution of n was examined, it was 


found in fact to be almost normal with a mean 
2 


Ss 
(n) of 72-3 and = small (= 0-005). The next 


step was an examination of observed distributions 
of x. One such distribution is given in Table 2, 


together with the appropriate calculated normal’ 


frequencies. It will be noted that agreement 
between observed and calculated frequencies is 
good, and acceptable according to the 2° test. 
This lends confidence to the formulae derived for 
the mean and variance of the expected distribu- 


tion, as well as for its shape. 


Table 2. Distribution of x observed frequencies after complete 
mixing, together with calculated normal frequencies as found 


from 
1 1 £\2] with 0-5 
= exp. | — 
oV2n 2 7 and a = 0:5/\/72°3. 
Observed Calculated 
x frequency frequency 
655-600 2 | 7 
620-655 2 2-7 12-1 
585-620 | | 
19 19-8 
31 82-2 
“480-515 32 37-9 
445—-480 32 31-1 
410-445 16 17-8 
9 75 
B40—375 3 li 2:3 10-1 
3°5--240 1 5 
491 
d. f. 6 
P {a4 > +91} 37 


4. INcomMPLETE MIXxING 
Theoretical considerations 
On the basis of the definition for complete mixing, 
it follows that mixing would be incomplete if the 
probability (@) of drawing red were not identical 
at all points in the mixer. Therefore, if mixing 
were incomplete and samples were withdrawn at 


points P,..... P,, then 2,(i=1....k) would 
be normally and independently distributed with 
—@,), 


means @, and _ variances where 


v 


v = E(n). The calculations are simplified if the 
variables 2, are transformed [8], [9] to new vari- 
ables z,; according to 


= 2 are sin 


For v fairly large, the effect of this is that at 
points P,..... P, variables z; would be observed 
which would be approximately normal with means 
1 
{, and variances -. 
Vv 

Some measure of the variation between the 
{, would be a measure of the degree of mixing, 
based upon the information obtained from 
points P,....P,. If mixing were vomplete 
= % (where is the trans- 
form of the total proportion of red in the batch) 

k 
and it follows that ¢=v (z; — would 
be a x variable. This could be used as a test 
criterion for the statistical hypothesis H, : 
C, = = as a test for 
completeness of mixing. 

If it were necessary to reject the H, hypothesis 
and therefore to conclude that mixing were not 
complete, the matter would have to be considered 
further. Although a “ mixing index” is useful 
in practice, it is quite arbitrary, and, as has been 
mentioned above, Lacey concluded that it was 
probably impossible to express the state of a 
mixture adequately by a single number. 

When mixing is incomplete, what appears to 
be important fundamentally is not the degree of 
mixing but rather the approach of the mixture 
towards the completely mixed state. As a means 
of determining this, ¢ may be used. 

A set of experiments was carried out to obtain 
data to test the theories outlined above. 


‘ 
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EXPERIMENTAL PROCEDURE 


The drum was filled as described previously, then 
attached to the stationary turntable at 30° to 
the horizontal. Definite sampling points were 


selected as shown in Fig. 4, the positions being 


HORIZONTAL PLANE 


located in relation to the horizontal fixed base of 


the mixer assembly. The batch was sampled 
before rotation was started, and the samples 
analysed by counting. The drum was then rotated 
for 20 seconds at 55 r.p.m., stopped and sampled. 
This mixing and sampling procedure was repeated 
At each 


sampling approximately 90 separate samples were 


a number of times as shown in Table 3. 


examined. The sand removed as samples in this 
test was not returned to the drum, as over the 
whole experiment it formed only 0-3°%, of the total 
batch. 


Table 3. Eaperimental mixing times. 


Sampling series No. Period of drum rotation — secs. 


Actual Cumulative 

I 0 0 

Il 20 20 
Ill 20 
20 60 
Vi 
Vil 150 
Vill 0 180 
IX tit) 2M) 
x 


Discussion oF ReEsuLts 

Adequacy of sampling 

In the section on the theory of incomplete mixing, 
it was assumed that & samples would be with- 
drawn at points P, Increasing the 
number of these sampling points would increase 
the amount of information about the state of 
the mixing, but for practical reasons this number 
must be kept small. The question, therefore, was 
whether a small number of points could give a 
true picture of the state of the mixing. 

It was considered reasonable to assume that 
after a period of mixing, the variation of @ in 
the mixer would be continuous. Further, there 
appeared to be no a priori reason for supposing 
that at any stage during the mixing process this 
continuous variation of @ would be a systematic 
fluctuation which would have the value @, at all 
the arbitrarily selected sampling points. Such a 
thing, while patently possible, could be regarded 
as extremely unlikely. Assuming that with the 
sampler used, the batch contained a total of N 
possible sampling points (N would be a very 
large number) and that at a proportion Q of 
these points @ = @, (where @, is the total propor- 
tion of red particles in the mixer), while at the 
others @ 
selected, the chance of selecting all these at points 
at which @ = @, is (Q)*. Now if S > 10, then 
even for Q as large as 0-5 this chance is < 0-001. 

In the experiment under consideration the 


é,, then if S sampling points were 


estimate of vy was 72-3 (cf. 3, Discussion of Results). 
After 60 seconds of mixing, samples were with- 
drawn from 89 points and calculations gave 
= 261-0. 
bution P { xf, 
critical level of 0-01, 17, had to be rejected. The 
calculation was repeated using only 10 of the 


According to tables of the 2* distri- 
> 261} < 001, so that, choosing a 


sampling points. These were selected so that at 
0 for 5 of the points and @ = 1 for 
As shown in Table 4, the result for 


66-7 and P { x}, > 66-7} 


zero time @ 
the other five. 
these 10 points was ¢ 
< 001, so that again H,, had to be rejected. 

Similar calculations after 150 seconds of mixing 
resulted in : 


84 points: ¢ = 79-7, P{ xf, > 79-7} > -60; 
accepted. 
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10 points: ¢ = 14-1, P{ xi) > 14-1} > -10; 
H,, accepted. 


In both cases the 10 points gave the same 
result as the much larger number of points. 
Although this is not a rigorous proof, it does give 
some confidence for using the smaller number of 
points and for regarding as reliable any con- 
clusions, concerning the degree of mixing of the 
batch, which may be based upon them. 


Process of Mixing 

Having decided that it was possible, on the basis 
of 10 sampling points to conclude whether or not 
mixing was complete at the time of sampling, 
values for ¢ were calculated for various times 
during the experiment. Results are given in 
Table 5, and have been plotted against time in 
Fig. 5. Information about the progress of mixing 
with time may be obtained from this plot by 
noting the rate at which the graph approaches 
k. Thus it appears that after 120 seconds, mixing 
was practically complete, as at that time ¢ 
approached its theoretical asymtotic value of 10.* 


Table 4. Typical calculation of ¢ (after 60 secs. mixing.) 


(z = 

| 2 are vin | 2 [(are sin 2) 

Point (are sin 
0,5,A 0-851 2-3400 | 0-7782 
| | 12694 | 0-2014 
2,4,B | 0-385 1-3387 | 0-2321 
3,1,A | 0-568 1-7072 | 0-1364 
4,.8,C | | 1-4263 | 0-1445 
5,6,B 0-468 15068 — 0-0640 
0,6,D | 0-250 12661 | — 03047 
17,A 0-377 1-3222 — 02486 
2,7,B | 0-422 14142 — 01566 
3,4,A | 0-493 1-5568 — 0-0140 

¢ = 72-3 {(-7782)? + (2014)? +...... + (-0140)?} =66-7 


** From tables [9]. + 2 are sin = 1-5708. 


* It can be shown that E(¢) =k + A 


where A=v Z(t; — &)*. 
i=1 
When all A=9, 


therefore when the system is completely mixed, E (¢) =k. 


Table 5. State of mixing after rotating the drum for various 


periods. 
2 
Time secs. ¢= — fo) 
i= 

0 1541-2 

20 311-7 

Ww) 124-0 

60 | 66-7 

90 36-1 

120 13-7 

159 14-1 

180 2-2 

240 7-7 

300 9-0 
1600;— 
wel 
soo | 
soof-| 


TIME 


Fig. 5. Progress of mixing with time based upon values 
of ¢ for 10 selected. 


5. CONCLUSIONS 


In considering the operation of solid-solid mixing, 
it was considered advisable to start with a highly 
simplified ideal case. Thus, the distinguishing 
feature of all the experiments discussed in this 
paper is that the particles to be mixed differed 
only in colour. Using this ideal case as reference, 
the theoretical concepts of complete and incom- 
plete mixing have been discussed and it has been 
concluded that the question of defining the 
completely mixed state is largely a statistical 
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one. It has been assumed that if, starting from a 
certain datum line, a set of mixing operations 
were performed and small samples were taken at 
various characterised points in the batch, and if 
this were repeated a number of times, always 
starting from the same datum line, then at each 
point, n, the number of particles in a sample and 
r, the number of red particles in the sample, could 
be represented by normal distributions. The 
proportion of red particles in the sample at each 
point may then also be represented by a normal 
distribution with a definite mean and variance. 
The definition of the completely mixed state 
follows as that which would obtain when the 
probability of drawing red at a certain point is 
the same at all points in the mixer. These basic 
assumptions have been confirmed experimentally 
and it has been shown, by using recognised 
statistical procedures, that a statistical hypothesis 
may be used as the test criterion for complete 
mixing. Furthermore, for application of this 
hypothesis, a comparatively limited number of 
samples is suflicient. In cases where the statistical 
criterion shows that mixing is incomplete, the 
progress of the mixture towards the completely 
mixed state may be readily followed. 


| 


It is clear that when “ ideal” conditions do 
not hold, the problem becomes much more com- 
plicated, and may necessitate certain simplifying 
assumptions. However, it is considered that the 
method of approach and the treatment used here, 
as well as the conclusions reached, will serve as a 
sound basis for further work. 
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NOTATION 


p(t) = Probability of an observation of size r. 
(") Binomial coefficient, "C,. 
~~ * Notation for moments about the mean and 
origin respectively [7] 
By. Bg = Coeflicients of skewness and kurtosis [7]. 
E (4) = * The expected mean value of x" [8]. 
P (a r Probability that 2 is greater than c. 


) 
2 Chi-square variable with k degrees of free- 
dom [8] 
var = Variance of 
d. f. Degrees of freedom. 
E (n). 
6 = Probability of drawing red at a certain point, 


Z; 2 are sin Va 
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The absorption of chlorine from air by solutions of olefins 
and iodine in carbon tetrachloride 


G. H. Roper 


School of Chemical Engineering, N.S.W. University of Technology, Broadway, Sydney, N.S.W. 
(Received 1 September 1953) 


Summary —The Steruens-Morais disc column is used tostudy the influence of chemical reaction 
upon the rates of absorption of chlorine from air by solutions of either cyclohexene or oleic acid 
with iodine in carbon tetrachloride. The increase in the liquid film coefficient due to chemical 
reaction between chlorine and cyclohexene was found to be independent of the concentration 
of cyclohexene, but no satisfactory correlation of the experimental results could be obtained. 

The fractional increase in the liquid film coefficient due to reaction between chlorine and 


oleic acid varied from zero to 9, and was found to be 
(ec + 0-00005)c 
$-1=% | OS 


where ¢ is the ratio of the liquid film coefficient to the coefficient in the absence of reaction, 
Cc is the concentration of the catalyst (iodine), ¢p is the concentration of oleic acid, and c; is the 
153 concentration of the dissolved chlorine at the gas-liquid interface. The units are lb. moles, feet 
and hours. The film theory is unable to predict the effect of the chemical reaction on the liquid 


film coeflicient. 


Résumé—Au moyen de la colonne A disques de l'auteur étudie influence 
d'une réaction chimique sur la vitesse d’absorption du chlore dans l'air au moyen de solutions 
diode dans CCl, contenant soit du cyclohexéne soit de l’acide oléique. 

L’accroissement du coefficient de transfert relatif au liquide accompagnant la réaction chimique 
entre le chlore et le cyclohexéne est indépendant de la concentration de ce dernier, sans qu'il ait 
été possible de trouver une corrélation satisfaisante entre les divers résultats expérimentaux. 

L’accroissement partiel du coefficient provenant de la réaction entre le chlore et l'acide oléique 
a varié entre 0 et 9 et peut étre représenté par : 

(ce . + 0-:00005)cp 


Cj 


¢—1=39 
ou ¢ est le rapport entre les deux valeurs du coefficient avec ou sans réaction, ce est la concen- 
tration du catalyseur (iode), cy, la concentration de acide oléique, C;, la concentration en Cl, 
dissous & l’interface gaz-liquide. Unités : Ib. moles, pieds, heures. 
La théorie du film de passage ne permet pas de prédire l'effet de la réaction chimique sur le 


coefficient de transfert relatif au liquide. 


INTRODUCTION another. StrepHeNs and Morris [6] designed a 

There has been little previous investigation on a | disc column which is suitable for small scale 

laboratory scale into continuous absorption of a | laboratory experiments and which gives for 

soluble gas followed by chemical reaction in the | physical absorption results similar to those 

liquid. Absorption in a wetted-wall column, | obtained in packed towers. 

although continuous, does not resemble absorption The dise column has been used by STEPHENS ~ 
in a commercial packed tower sufficiently closely | and Morris to study the absorption of chlorine é 


to enable reliable predictions to be made of the 
mass transfer coefficients in the commercial 
tower. In a packed absorber the liquor is mixed 
each time it flows from one piece of packing to 


by solutions of ferrous and ferric chlorides. With 
a total iron concentration equivalent to 28-7 lb. 
of FeCl, per cu. ft., the ratio of the ferrous 
chloride concentration (c,) to the concentration 
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dissolved, but unreacted, chlorine at the gas- 
liquid interface (c,) was varied from zero to 
700. The fractional increase in the liquid film 
coefficient due to chemical reaction was found 


to be 
083 
0-75 


The absorption of chlorine from air by solutions 


(1) 


2-ethyl hexene-1 and iodine in carbon tetra- 
chloride has been studied using a similar dise 
column (3). The ratio of the concentration of the 
olefin (c,) to the concentration of dissolved 
chlorine at the interface (c,;) was varied from zero 
to 300. The increase in the liquid film coeflicient 
due to chemical reaction was found to be 


¢—1 | 


os 


(2) 


be, 


where b varied linearly with the concentration of 
the catalyst (iodine) and was taken to be pro- 
portional to the specific reaction rate for the 
addition of chlorine to 2-ethyl hexene-1 under the 
conditions of the absorption. In the absence of 
iodine, ) was 0-08 and was 1-68 when the iodine 
concentration was 0-001 Ib. moles per cu. ft. 

The kinetics of the addition, in carbon tetra- 
chloride, of chlorine to 2-ethyl hexene-1, to oleic 
acid, and to cyclohexene have been investigated 
[4]. The of the addition to 2-ethyl 
hexene-1 and to oleic acid may be expressed 


(3) 


dt 


where c, and c, refer to the concentration of 


the chlorine and olefin, respectively. The value 
of the specific reaction rate k. at 70°F is 155 
l. (g. mol.) (see.) or 9 10° cu. ft./(lb. mol.) 
(hr.). The kinetics of the addition to cyclohexene 


may be expressed 
4 (4) 


the reaction being second order with respect to 
the chlorine concentration and zero order with 
respect to the cyclohexene concentration when 
there is an excess of cyclohexene in the liquor. 
The value of the specific reaction rate, at 70°F, is 


The absorption of chlorine from air 


640 |. /(g. mol.) (sec.) or 87 x 10° cu. ft./(Ib. mol.) 
(hr.). 
EXPERIMENTAL PROCEDURE 
The apparatus and procedure have been pre- 
viously described fully [8], [6]. The absorption 
unit is a StrepHens-Morris disc column whose 
major dimensions are : 


Number of discs 81 
Disc diameter (mean) 1-475 cm 
Dise thickness (mean) 0-45 cm 


2.50 cm 
0-122 ft. 
0-052 ft. 
89% 
0.184 ft. 


Tube diameter 

Mean perimeter for liquid flow 
Equivalent diameter for gas flow 
Free space (dry) 

Absorption surface (dry) 


The dises are threaded edgeways on a 2 mm 
glass rod and successive discs are maintained at 
right-angles. 

The overall coeflicients were calculated from 
the equation 

kK (5) 


and were corrected for the resistance in the gas 


N ,/AAp 


G 


phase in order to evaluate the liquid film coeffi- 


cients. The gas film coeflicients were calculated 
from 


and the liquid film coefficients in the absence of 
a reaction from 
Hk, 


0.0041 (7) 


For each olefin (i.e., oleic acid or cyclohexene) 
the first series of runs were made without using 
the iodine catalyst. The concentration of the 
physically dissolved chlorine in the liquor leaving 
the discs was determined from the concentration 
of chlorine in a hydrochloric acid solution placed 
in the liquor run-off tube. The total concentration 
of chlorine in the liquor was determined by the 
decrease in the olefin concentration after sufficient 
time had been allowed for all the chlorine to 
react. There was no smell of hydrogen chloride 
from the liquor, indicating that there was no 
substitution. 
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Table 1. Absorption by solutions of cycloherene and iodine in carbon tetrachloride. Constant air rate = 0-790 c.f.m. 
| 
Liquor 
— -- - - Gas | 
| 108 cg | 108 | 10% cc | mean 
(ft.) in | out free | total | | 
1 79-1 82 59 0-32 0-65 2-52 | 0-57 0 | O44 | 0-063 0-104 0-131 
2 122 68 0-89 2-53 10 | 0-93 | 0-146 0-129 0-171 
3 845 | 70 | 52 | 0-80 | 0-51 69 | 0-49 0 040 | 0-051 | 0-101 | 0-127 
4 78-0 71 54 O42 | 1:13 66 | 10 | 0 0-84 | 0-108 | 0-084 | 0-103 
5 70-6 65 53 | 0-31 1-59 10:3 | 12 | 0 101 | 0-137 0-081 0-098 
6 733 | 66 | 52 | O36 | O72 | 106 | O71 | O 0-68 | 0065 | 0-074 | 0-089 
7 68-5 | 73 | 66 201 | 31 | O78 0-65 0-77 | 0-168 | 0-124 | 0-168 
s | 431 | 73 | 55 1-46 32 | 050 | 065 | 0-47 | 0076 | 0-092 | 0-116 
9 | 118 73 56 0-65 3-9 0-41 0-65 040 0-094 | 0-137 | 0-190 
10 | 185 «|| 73 63 160 | 31 | 125 | 0-65 115 | 0-293 | 0-141 | 0-183 
| | | 


With iodine in the liquor, the presence of the 
interhalogen compound prevented the use of the 
method for determining the concentration of the 
physically dissolved chlorine in the liquor leaving 
the discs. The ratio of the physically dissolved 
chlorine to the total amount of chlorine absorbed 
was estimated from a graph in which this ratio 
was plotted as a function of N,/cJ° in the 
absence of a catalyst. 


Resu_ts anp Discussion 
Absorption by solutions of cyclohexene 


The experimental data for this absorption are 
given in Table 1. Six runs were made without 
catalyst at liquor rates from 70-6 to 122 Ib. /(hr.) 
(ft.), and four runs were made at liquor rates 
from 43-1 to 185 lb./(hr.) (ft.) with an iodine 
concentration of 0-65 10° Ib. moles per 
cu. ft. 

The addition of a small amount of cyclohexene 
to the liquor passing over the discs increased the 
liquid film coefficient, but further additions did 
not increase the liquid film coefficient (Runs 3, 
6). This agrees with the kinetics expressed by 
equation (4) in which the rate of disappearance 
of chlorine is independent of the concentration of 
cyclohexene provided that there is sufficient 
olefin present to react with all the chlorine. The 
catalysed absorption has the same characteristics 
which suggests that the rate of the catalysed 


A 


reaction is similarly independent of the con- 
centration of cyclohexene. 

No satisfactory correlation of the data could 
be obtained. Based on either the film theory or 
the penetration theory, an elementary treatment 
of absorption followed by second-order disappear- 
ance of the gaseous solute indicates that the 
liquid filn: coefficient should increase with in- 
creasing values of c;. The data of Table 1 (runs 1, 
8, 4, 5, 6 and runs 9, 10) show that c, has the 
opposite effect on the coefficient. Comparable 
results have been obtained by Preaceman [2, 
p. 118] for the desorption of chlorine from water 
by air. He found that the film theory fails badly 
in explaining his experimental results. The most 
significant aspect of the failure is that the film 
theory predicts that the liquid film coefficient 
decreases with increasing concentration of 
chlorine, while it was found experimentally that 
the liquid film coefficient remains constant for 
low concentrations and then increases sharply 
with concentration. 


Absorption by solutions of oleic acid 


The experimental data for this absorption are 
given in Table 2. Ten runs were made without 
using catalyst at liquor rates from 50-2 to 241 
Ib. /(hr.) (ft.) and seven runs were made at liquor 
rates from 54-5 to 118 lb./(hr.) (ft.) and with 
iodine concentrations from 0-09 x 10° to 0-96 
x 10-8 Ib. moles per cu. ft. 
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Table 2. Absorption by solutions of oleic acid and iodine in carbon tetrachloride. Constant air rate = 0-790 c.f.m. 


Liquor | | 

Run | Rate | Temp. °F. | 108 | % Cl, 10° Ny 
No. /(hr.)|—— 10% | 16% | 16° e¢ Kg Hk, 

(ft.) in out free | total | | out in liquor | gas 

| 

1 83-6 | 65 56 | O87 | 32 214 | 22 | 0 | 19 28 | O83 | | 0-092 | O-119 
2 60-6 | 62 59 | O14 | 47 «| «16-7 15 | 0 | 13 16 | O85 | OF | O134 | 0-182 
3 53-3 | 63 54 | O11 | 21 180 | O8 | O 06 OF | O13 | OL | O115 | 0-154 
4 61-1 | 64 62 | 0-27 | 50 167 | 29 | Oo | 28 27 | 0-87 05 | 0-086 | 0-106 
5 | 241 71 68 | 0-51 | 38 185 | 17 | O 22 B81 | 118 1-1 | | 0-42 
6 | 218 70 70 | 069 | 50 198 | 31 | o | 31 45 | 1-33 (16 | 0208 | 0-81 
7 | 202 73 | 65 | 0-25 | 19 | 22 |) 14 | O 10 «14 | O47 | O5 | 0-228 | O38 
8 714 | 69 56 | O13 | 18 88 | 09 0 O7 O8 | O16 | OL | O112 | 0-148 
” 60-4 | 64 56 | O18 | 40 | 209 | 12 ft) 10 12 | O29 | O2 | 0138 | 0-200 
10 50-2 | 64 64 | 094 | 90 | 187 | 46 0 04 O48 | O55 0-083 | 0-102 
11 | 118 69 62 3-2 26 | 15 oss | 14) 619 | 046 | 06 | 0-165 | 0-25 
12 715 | 69 | 51 10 40 | O89 | O22 | O4 OF | 0-09 | 0-146 | 0-21 
13 65-6 | 69 56 24 34 | O76 | 045 | OF O8 | O19 | | O41 0-20 
14 616 | 68 54 19 112 | O53 | O22 | OF O58 | O14 | 0-156 | 0-23 
15 545 | 69 B44 | 028 | 153 | O24 | 095 | 04 0-4 0-019 0-27 0-68 
16 80-2 | 72 63 22 | 143 | 15 009 | 10 13 | O22 0-152 | 0-22 
17 902 | 70 | 64 50 | 6-4 0-9 096 | 12 15 | O55 0-23 0-44 


The film theory of gas absorption has been | gaseous solute with one component of the liquor 
examined by PeaceMAN [2, pp. 206-31]. The effect | is determined by two dimensionless groups : M 
of an irreversible second-order reaction of the | or k,c,D,/(k,')*; and q or Dyc,/Dgey. The 


Table 3. Absorption of chlorine by solutions of oleic acid. 


| 
| ¢-1 
Run r | | | Hky/ VM q theory 
No. ezpt. | 
| film pen 
1 83-6 21 2-2 0 0-119 0-091 1-7 2-3 0-3 0-5 
2 60-6 16-7 15 0 0-182 0-073 66 | 0-20 15 3 5 
3 53-3 18-0 0-8 0 | 0-154 0-066 70 0-099 13 4+ 
4 61-1 16-7 2-9 0 | 0-106 0-073 66 | 0-39 0-5 2 
5 | 261 18-5 1-7 0 O42 0-191 24 | 0-20 12 1-4 
8 714 8-8 0-9 0 | ©1438 0-081 40 | O23 0-8 4 
9 60-4 20-9 12 0 O20 0-076 66 | 0-13 16 4 
10 50-2 18-7 46 0 0-102 0-064 2-4 | O55 0-6 2 
11 118 2-6 15 | 045 0-25 0-116 4-6 13 12 | 0-8 
12 71-5 40 0-39 0-22 0-21 0-080 6-2 0-22 | 16 | 3 3 
13 65-6 84 | O76 | 0-45 0-20 0-077 81 | O50 | 16 2 
14 61-6 112 0-53 022 | 0-23 0073 | Ol | 22 6 9 
15 54-5 15-3 024 | 096 | 0-68 0-068 28 0-035 9 17 | 
16 802 | 143 15 009 | 0-22 0-088 78 | 023 15 3 | 4 
17 90-2 6-4 0-9 096 | O44 0-096 13 0-31 3-6 3 | 3 
| 
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theoretical treatment of PeEacEMAN is limited by 
the assumption that none of the gaesous solute 
is unreacted in the main body of the liquor. This 
requirement is approximated in those runs in 
which the concentration of the free (unreacted) 
chlorine in the liquor leaving the discs is small 
compared to the total concentration of absorbed 
chlorine. In Table 8 there are given the recalcu- 
lated data from the runs using oleic acid in which 
the ratio ¢ (= k,/k,') is predicted by the theory 
to be greater than 3; the fraction of chlorine 
unreacted is expected to be small for these runs. 
The results previously published using 2-ethyl 
hexene-1 (8) have been similarly re-calculated and 
are given in Table 4. The diffusivities of the 
solutes, for use in calculating the dimensionless 
groups, were estimated by the method of ARNOLD 


| 


[1]. The calculated values of the fractional increase 
in the liquid film coefficient (¢ — 1) were read 
from the graphical solution of Peaceman. 
Perry and Picrorp [7] have published a 
solution for unsteady-state absorption and second- 
order chemical reaction (penetration theory). By 
assuming equal diffusivities, the effect of the 
chemical reaction is determined by two groups, a 
dimensionless time k,c,t and c,/c;. For certain 
limiting cases, the effect of unequal diffusivities 
is given by replacing by D,c,/D ec; In 
evaluating ¢ from the graphs of Perry and 
PicrorD, it was assumed that the effect of D, 
and D, is the same as in the limiting cases. 
PEACEMAN [p. 197] has shown that the time, t, is 


t = 4D ,/3-14 (k,')? 


Table 4. Absorption of chlorine by solutions of 2-ethyl hexene-1 
| 
Run 10° cp VM 108 ¢; q 108 | theory 
No. expt. 
| | film | pen. 
2 47 2-5 49 19 ons | 0 | | os 
3 25 46 0125 | | os | 06 
4 | 23 | | os | 0-8 
5 31 64 | 0-96 0054 | 0-26 0 | 2 
6 30 6-3 | 22 0-124 | O147 | O | O6 4 2 
8 42 (58 3-2 0158 | | 0-8 
9 20 0-77 0-188 020 | #10 | | 
10 58 (62 11 0-33 0-25 0-20 16 | 
6-8 0-45 0-11 033 0-20 25 4 | of 
12 9-1 0-19 0-037 0-88 | o50 | 77 8 3 
13 6-2 o-4 15 O41 =| 028 | 0-50 20 3 
14 5-0 4 2-3 0-80 0-22 | 0-50 13 1-2 
15 13 19 0-53 0-073 0-61 1-04 5:5 11 14 
16 34 31 0-12 0-0062 2-0 | cl 29 
17 82 O87 0-020 13 104 | 22 
18 2 1-2 0-073 0-61 104 54 12 
19 27 28 2-2 0-14 0-53 | 1-04 45 7 7 
20 21 12 0-10 | 0-081 | 7 
21 18 W 1-2 | 0-12 oe | o7 | 46 5 2 
22 12 6-9 0-75 O11 067 | 0-42 31 | 4 9 
23 66 14 | O22 0355 | O42 24 | 5 
24 nu 24 | 087 | 0 06 2 
26 13 2-4 036 | 0048 | Ott | 0 1-9 2 
27 12 6-7 24 | 0-36 0-46 0-42 18 3 3 
28 9-5 8-0 2:5 | 0-46 0-56 | 0-77 24 2 2 
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and the other dimensionless group is thus 
4k,c,D ,/8-14 (k,')*, or 4M /3-14. 

For both the film theory and the penetration 
theory, the calculated values of (¢ — 1) are of 
the correct order. Both theories fail in that they 
predict that M is an important dimensionless 
group. If M were important, the correlations of 
the previous paper [3, Figs. 5, 6 and 7] would 
depend on the liquid film coefficient, k, , and hence 
on the liquor rate. Similarly, the points plotted 
in Fig. 1 of this paper would be scattered accord- 
ing to liquor rate. The correlation of SrEPHENS 
and Morxis for the absorption of chlorine by 
ferrous chloride solutions similarly shows that 
M is not important as their correlation of (¢ — 1) 


against c,/c; is independent of liquor rate.* 


The reaction kinetics for the addition of chlorine 
to oleic acid in carbon tetrachloride are, within 
experimental error, the same as the kinetics for 
the addition to 2-ethyl hexene-1. Equation (2) 
for the absorption of chlorine by solutions of 
2-ethy! hexene-1 and iodine in carbon tetrachloride 
may be re-written 


(8) 


¢—1=39 


In Fig. 1, the data from Table 1 has been 
plotted according to equation (8). That the data 
for absorption by solutions of oleic acid is in 
excellent agreement with the data for absorption 
by solutions of 2-ethyl hexene-1 is shown by the 
grouping of the experimental points about the 
solid line representing equation (8). 


O tc =0 

7 ¢¢ = 0-09 x 10° Ib. mol /ft.* 
x = 022 x 10° Ib. mol /ft.* 
+ ¢¢ = 045 x 10° Ib. mol /ft.* 
tc = 0-96 x 1078 lb. mol /ft.* 


[> 


* “In a recent paper, VAN Kreveten and Horryzer 
[8] have reported a graphical design method for gas-liquid 
reactors. Their method is based on the film theory, but is 
inapplicable to a disc column as the group M is used. The 
dise column more closely resembles a commercial packed 
tower than does either of their experimental units, and 
hence their design method should be used with considerable 
caution for commercial towers.” 
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Fig. 1. The fractional increase in the liquid film coefficient 
due to the chemical reaction ofchlorine in solutions of oleic 
acid and iodine in carbon tetrachloride. 


CONCLUSION 


The film theory has been previously criticized 
because it fails to predict the magnitude of the 
absorption coefficients in the absence of reaction, 
because it predicts that these coefficients are pro- 
portional to the diffusivity of the gaseous solute, 
and because it fails to explain the effect of con- 
centration on the rates of desorption of chlorine 
from water [2], [5]. 

The film theory fails to explain the effect of the 
olefin concentration on the rates of absorption of 
chlorine by solutions of cyclohexene and iodine in 
carbon tetrachloride. The theory predicts that 
the fractional increase in the liquid film coefficient 
for the absorption of chlorine by solutions of 
2-ethyl hexene-1 or oleic acid with iodine in carbon 
tetrachloride is affected by the liquor rate, whereas 
the experimental values of the fractional increase 
are independent of the liquor rate. 
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NOTATION 

A = surface area for mass transfer, sq. ft. 

b = a parameter defined by equation (2). 

¢4 = concentration of chlorine, lb. moles per cu, ft. 

Cp = log-mean concentration of the liquid phase re- 
actant, lb. moles per cu. ft. 

‘c = concentration of catalyst (iodine), Ib. moles per 
cu, ft. 

c = concentration of dissolved, but unreacted, chlorine 


at the gas-liquid interface, lb. moles per cu. ft. 
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D4 = molecular diffusivity of chlorine, sq. ft. per hr. 
Dg = molecular diffusivity of liquid phase reactant, sq. 


ft. per hr. 

Gy = mass velocity of the gas stream, Ib. moles /(hr.) 
(sq. ft.). 

H = Henry's law coefficient, Ib. moles /(cu. ft.) (atm.). 

kg = specific reaction rate defined by equation (3), cu. 
ft./(Ib. mole) (hr.). 

hy = specific reaction rate defined by equation (4), cu. 
ft. /(Ib. mole) (br.). 

kg == gas film coefficient of mass transfer, Ib. moles /(hr.) 
(sq. ft.) (atm.). 

ky; = liquid film coefficient of mass transfer, ft. per hr. 


ky’ = liquid film 2oefficient of mass transfer in the 
absence of reaction, ft. per hr. 


Kg = overall coefficient of mass transfer, Ib. moles /(hr.) 


(sq. ft.) (atm.). 

M a dimensionless group, k, ¢g D4/(k,')?. 

N4 = the rate of absorption of chlorine, lb. moles per 
hr. 

Ppm = the mean partial pressure of air in the gas film, 
atm, 

Ap = the logarithmic mean between the driving forces 
at each end of the column, atm. 

q = a dimensionless group, D 4 ¢;/Dz cp. 


Reg = the Reynolds’ number for flow in the gas stream. 
Scqg = the Schmidt number for the gas phase. 

t = time hr. 

r = the rate of flow of the liquid layer, Ib. /(ft.) (hr.). 
the ratio of the liquid film coefficients = k, /k;’. 
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A further study of the hydraulic cyclone 


D. F. 
Atomic Engergy Research Establishment, Harwell, Didcot, Berks. 


( Received 30 October 1953) 


Summary—By a technique involving injection of closely sized batches of Perspex spheres 
into the feed of a 3 in. hydraulic cyclone, the effect of several variables on solid elimination 
efficiency have been investigated. 

A } in. feed diameter results in maximum efficiencies on this cyclone for a wide range of 
conditions. It is considered impossible to derive simple power relationships to describe the 
effect on particle elimination efficiency of change in feed diameter, overflow diameter and through- 
put, to cover a range including the optimum conditions. 

The importance of turbulent mixing in the feed section due to shock effects and of the short 
circuit flow down the outside wall of the vortex finder have been demonstrated. 


Résumé—L auteur étudie influence de plusieurs variables sur lefficacité de la séparation des 
poussiéres en alimentant un cyclone hydraulique de 3 pouces avec des charges constituées de 
sphéres en Perspex (plexiglas) dont les diamétres restent compris entre des limites étroites. 

Pour lalimentation de ce cyclone, un diamétre de } de pouce donne une efficacité maximum 
dans une échelle étendue de conditions. Pour un cyclone donné, il parait impossible de prédire 
lefficacité du classement par une relation générale faisant intervenir, avec des exposants 
simples, chacun des facteurs : diamétre de l'alimentation et du trop-plein, débit, etc. . . pour 
embrasser les divers modes de fonctionnement, de par et d'autre des valeurs les plus favorables. 


Ce travail souligne l'importance du mélange turbulent dans la section d’alimentation et 
du court-circuitage vers le bas de la paroi externe du conduit central. 


1. INTRODUCTION 


The extensive literature on the hydraulic cyclone 
contains little information about carefully con- 
trolled experiments designed to give an accurate 
assessment of the effect of the several cyclone 
variables on particle separation (or elimination) 
efficiencies. Most of the work reported involved 
the use of heterogeneous particles of various shapes 
and was carried out with a mixture of particle 
sizes under such conditions that an indeterminate 
amount of interference of the settling of particles 
in any particular size range could be exerted on 
particles of another size range. 

Several methods of comparing efficiencies have 
been used. Dantstrom [1] uses d,,, the particle 
size in microns which reports 50% by weight to 
overflow and underflow. Fitcn and Jounson [2] 
prefer to use a particle size which reports almost 
completely to the underflow (say 95°). A more 
reliable method involves the construction of 
efliciency curves which give the efficiency of 
elimination of the whole range of particle sizes 
present in the cyclone feed. 


One aim of the present work was to assess the 
relative merits of these methods of summarising 
cyclone efficiencies. Another was to examine the 
importance of short circuiting of particles to the 
overflow, a phenomenon noted by the present 
author in a previous paper [8]. 

In an attempt to decrease the number of 
variables, pulp densities used in the present study 
were kept very low and throughout the testwork 
dilutions were effectively infinite. To eliminate 
the effect of shape factors and differing specific 
gravities of particles, almost perfect Perspex 
spheres of constant specific gravity were used. 


2. Expertmentat Detai.s 
Apparatus 

Fig. 29a shows the cyclone to full scale. Feed 
openings and underflow diameters were varied 
by inserting interchangeable cyclone sections ; 
vortex finders, with jin. thick walls, were 
available such that lengths and overflow diameters 
were variable. 


Fig. 29b provides a diagram of the apparatus. 


254 


VOL. 
2 
1953 
| 
4 
6 


D. F. Kexsauz : A further study of the hydraulic cyclone 


For low pressure experiments the conical feed 
tank, fitted with overflow, was raised to the 
required height to give a constant static head. 
For pressures beyond the range of the head 
available the conical feed tank was connected to 
a centrifugal pump (or two in series) and the 
feed pressure adjusted by the setting of the valve 
between the pump and the pressure gauge. 


Preparation of sized fractions of Perspex spheres 
A miscellaneous mixture of Perspex spheres was 
closely sized by repeated wet screening in the 
presence of Aerosol OT as wetting agent. Tyler 
44/2 series screens were used down to 825 mesh. 
The minus 825 mesh fraction was further split 
into sized fractions using conventional elutriation 
in a multicone laboratory classifier. The classified 
fractions were improved by repeated sedimenta- 
sion in beakers in the presence of Aerosol OT until 
tatisfactory fractions were obtained. 

Representative samples of all fractions were 
examined under a microscope fitted with a 
graduated eyepiece, several hundred grains 
counted and their diameters measured. The 
effective average diameter of each fraction was 

3 
calculated as J “< and for all graphs shown 
in this report these average diameters were used 
to describe the particular fraction used for a given 
test. 

Specific gravity determinations of samples of 
the various sized fractions showed that the specific 
gravity was 1-18 in all cases. Inclusions of air 
were very rare and a negligible fraction of the 
particles were non-spherical. 


Test Procedure 

The desired combination of cyclone variables was 
obtained by selection of the appropriate feed 
opening, vortex finder and underflow extension 
cone. The head, or pump delivery, was adjusted 
to give the correct operating pressure and under- 
flow and overflow water rates were measured. 

A sample of the particular batch of Perspex 
spheres, always stored under water in the presence 
of Aerosol OT to give complete wetting, was 
obtained by pipetting off a selected volume from 
the main sample, which had been well mixed by 


agitation, and placed in a glass separating funnel 
fitted with stopcock and long delivery tube. After 
agitation the stopcock was opened, the delivery 
tube filled with slurry and closed at the end with 
one finger. The stopcock was then closed. 

The end of the delivery tube was immersed in 
the water in the feed cone, finger removed quickly 
and the tube inserted so that the open end was 
located in the throat of the feed funnel. The 
stopcock was then opened and all solids washed 
down with water containing Aerosol OT. These 
precautions were necessary to ensure that no air 
bubbles were introduced into the system. 

After a time interval, found by a preliminary 
run with no solids but using a dye, a double 
deflector was moved so that the underflow and 
overflow streams were collected in separate con- 
tainers before any solids appeared in the trans- 
parent cyclone. Collection was continued until 
no solids were visible within the cyclone and the 
deflector returned to its initial position. 

The Perspex spheres in each container were 
allowed to settle under quiescent conditions for a 
time known to be sufficient to permit settling of 
all particles in the particular fraction used. The 
supernatant water was carefully syphoned off. The 
solids were washed with a minimum of water, 
containing wetting agent, into smaller glass 
containers and the procedure repeated until the 
underflow and overflow solids were obtained as 
samples containing approximately 50°, water by 
weight in small beakers of known weight. 

The final stage involved oven drying to constant 
weight at 80°C, followed by accurate weighing 
to give the weights of the dry solids. 

The total weight of solids used in any test 
never exceeded 5 or 6 gm and the total volume 
of water collected often exceeded 10 litres. 


8. GRAPHICAL PRESENTATION OF RESULTS 


After calculation of the percentage of the total 
weight of the spheres which reported to the 
underflow, graphs of percentage eliminated to 
underflow against the average particle size, 
measured previously, were drawn. The point 
corresponding to zero particle size was that given 
by the percentage of the total water volume flow 
delivered to the underflow. 
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A knowledge of cyclone operation leads to the 
conclusion that if 2°, is the fraction of the total 
water which is discharged through the underflow 
then, independent of any centrifugal forces acting 
on particles, r°%, of all particles must leave through 
the underflow. The additional percentage of 
particles of a given size which are discharged 
through the underflow is a measure of the 
efficiency of elimination of solids due to the 
cyclone action. 

If y% is the experimentally determined fraction 
of particles of a given size eliminated via the 


underflow, then (y — x) % is the fraction elimi- 
nated from a total of (100 — z) % due to cyclone 


action. Consequently eas x 100 is the 


corrected percentage of particles of a given size 
eliminated through the underflow as a result of 
centrifugal forces within the cyclone. 

The graphs provided in this report show the 
relationship between this corrected percentage 
eliminated via the underflow and particle size. 

Not all uncorrected curves are presented but, 
if necessary, they can be constructed by back 


Table 1 
Conditions for Series 1-6 

P.D. O.P.D. | UPD. P T.F. U.P. Vol. 
(ins.) (ins.) (ins.) | (ins.) lbs /sq. in. litres /min. litres /min. % ULF. 

i | i 4 i 39 23-77 3-23 13-6 
| 18 16-80 2-47 
Series 1 | 10 12-97 2-00 15-4 
| 6 10-47 1-74 16-6 
; } 4 | i 20 13-83 1-57 11-3 
| 12 11-10 1-39 12-5 
8 9-34 1-29 13-8 
Series 2 6 8-23 1-20 14-6 
| 4 6-95 1-09 15-7 
| 2 5-19 0-95 18°3 
A 2 12-31 3-30 26-8 
| 20 9:34 2-63 28-2 
Series 3 | | 12 7-54 2-21 29-3 
| | | 8 6-39 1-93 30-2 
| ; | 2 | 15-82 1-54 
20 11-80 1-32 11-2 
Series 4 | | 12 957 | 118 12-3 
8-03 1-09 13-5 
6 702 1-02 145 
i 2 18-24 0-55 3-1 
20 13-81 0-53 3-9 
Series 5 12 11-01 0-53 48 
| 9-21 0-51 5-6 
i 2 20-42 0-21 1” 
20 15-32 0-28 18 
Series 6 12 11-87 0-27 2-2 
9-86 0-27 27 


F.D. = Feed Diameter. O.F.D. = Overflow Diameter. 


U.F.D. = Underflow Diameter. V.F.L. = Vortex Finder 


length. P. = Pressure. T.F. = Total Flow. U.F. = Underflow. 
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calculation since in all cases the percentage water 
flow to the underflow is provided. 


4. Derarts or Testwork 


Due to the large numbers of experiments which 
would be necessary to undertake a complete 
survey of cyclone operation, selected combinations 
of the numerous variables were chosen in order to 
throw some light on the effect of the more 
important ones on particle elimination efficiency. 

During most of the work the temperature of the 
water used during the experiments was not 
recorded, as the available information indicated 
that the mains water temperature was practically 
constant throughout the year within one or two 
degrees centigrade. However, experiments origin- 
ally carried out during the winter months did not 
give the same results when repeated in the 
summer. In all cases the deviations could be 
accounted for by the increase in viscosity of water 
with decrease in temperature. 

Most of the experiments within any one of the 
above sections were carried out during conditions 
when deviations of temperature were slight but 
precise comparison of the results of one section 
with those of another is difficult due to ambient 
temperature differences reaching approximately 

With very rare exceptions duplication of 
experimental results under exactly the same 
conditions was excellent. Consequently, apart 
from preliminary runs, no further tests were run 
in duplicate. 


The Effect of Change in Feed Pressure 
Series 1 to 6 


For initial experiments 4 in. long vortex finders 
were used but, after work had been completed 
which showed the effect of change in vortex 
finder length, it was decided to use vortex finders 
of 2 in. length. Consequently the first two series 
reported using gin. and gin. diameter feed 
openings involved 4 in. long vortex finders, } in. 
overflow and 4 in. underflow diameters. For the 
next four series 2in. long vortex finders were 
used and in all four cases the feed diameter was 
tin. but overflow diameters were ¥,, }, # and 
# inches respectively. 


Table 1 gives a summary of the conditions and 
water volume flows for the six series. 

Figs. 1 to 6 show graphs of corrected weight 
percent to underflow against particle size. 
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Fig. 2. Effect of feed pressure on particle elimination : 
Results of Series 2. 
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Fig. 3. Effect of feed pressure on particle climination : | pj, 


5. Effect of feed pressure on particle elimination : 
Results of Series 3. 


Results of Series 5. 
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Fig. 7 gives the relationship between d,, and 
feed pressure, and Fig. 8 between total through- 
put (G) and feed pressure. (In constructing Fig. 7 
small corrections were made for changes in water 
viscosity with temperature. On the assumption 
that Stokes’ Law governs the settling rates, 
dy where = viscosity.) 
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2 
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microns 


Fig.7. Relationship between feed pressures and d,, values: 
Series 1-6. 


Fig. 7 shows that for series 1, 3, 4, 5, and 6 


P «x ~, and for series 2 
1 
P 
50 


From Fig. 8 it can be seen that for all series 
G a or P x G**, 
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Fig. 8. Relationship between feed pressure and total 
throughput : Series 1-6. 


Hence for series, 1, 8, 4, 5 and 6 
d 1 and for series 2 
ep OF Gre and for series 2 


1 

50 
Further experiments are needed to throw more 
light on this phenomenon but it may be concluded 
that for any given combination of cyclone 
dimensions there exists a simple power relation- 
ship between feed pressure, or throughput, and 
deo, 
power is not constant. 


and that the limited results indicate that the 


The Effect of Change in Feed Diameter and 
Shape of Feed Opening 


The Effect of Change in Feed Diameter 
For the first series feed pressure was maintained 
constant at 6 lbs. /sq. in., overflow diameter } in., 
vortex finder length 4 in. and underflow diameter 
} in. Selected particle sizes were passed through 
the cyclone for a range of feed diameters. 
The second series was a repetition of the first 
but with a constant feed pressure of 15 Ibs./sq. in. 
For the third series the length of the vortex 
finder was decreased to 2 in. but all other variables 
were maintained as for series 2. 


59 


j +—— + + —+ — 
5 | + 
7 | 
| 
10 —+— + | 
| 
| 
‘ 
IL. 
) 
53 
| 
| 
| 
| 
: 


D. F. Kersati : A further study of the hydraulic cyclone 


The fourth series involved two sets of condi- a ere. 


tions selected at random to give wide changes in 
overflow and underflow diameters. 

The fifth series involved a repetition of series 2 
and 8 but using vortex finders j in. in length. A 
further test was made in which no vortex finder 0 


i Te) 


was used, the overflow being a circular hole in | 
the top of the cyclone (i.e. zero length vortex 60 

finder). 
Figs. 9, 10, 11, 12 and 18 summarize the results : so 

of the respective series. 5 
z 40r 

ee = 


PARTICLE 
Size 


> microns 
wo] | 
/ \ v 2.0 


| T 
5 / FEED DIAMETER 
2 / 
ze : Fig. 10. Effect of change in feed diameter : second series. 
3 
el*lele iz 
>| > imucrons 
| \ \ 
T 303 \ 
4 i i 
00! 
Fig. 9. Effect of change in feed diameter : first series. \ 
\ 
The weight percents graphed are those obtained A 
by experiment and are uncorrected. No significant # . 
alteration of the shape of the curves was obtained : 
using corrected results. Consequently, for con- 3 3] 
venience, actual results were used. 2 | LJ’ |mcrons 
Inspection of Figs. 9 to 13 reveals that, without | 
exception, for a wide range of conditions, an a | at goog 
optimum feed diameter exists with an approximate & 1 
Figs. 9, 10 and 11 show the marked increase e466 5 
in efficiency for all particle sizes as the feed 
diameter decreases from § in. to } in. Fig. 11. Effect of change in feed diameter : third series. . 
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Fig. 12, the results of conditions selected at 
random, shows that within the accuracy of the 
experiments, a feed opening of } in. diameter 
gives optimum particle elimination for overflow 
diameters of }in. and } in. 
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Fig. 12. Effect of change in feed diameter : fourth series: 


No further experiments were undertaken but 
it is inferred that, for all practical purposes when 
using a 8 in. cyclone, optimum elimination of all 
particles is obtained using a feed diameter of 
din. to in. 

Additional conclusions may be drawn from a 
consideration of the curves for 49 microns particles 
in Figs. 10, 11 and 13, which throw considerable 
light on the effect of vortex finder length on the 
efficiency of elimination of relatively coarse 
particles. However, this aspect is dealt with later. 


The Effect of Change in Feed Diameter on the Rate 
of Injection of Momentum 


In view of the interesting results given in the pre- 
vious section a detailed study was made of total 
volume flows (through-puts, G) for a wide range 
of cyclone conditions. The procedure involved 
selection of a given feed pressure, overflow 


diameter and underflow diameter and measure- 
ment of total flows for a range of feed diameters. 
The vortex finder length was maintained constant 
at 4 in. 

Average feed velocities (V,,) were calculated 
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Fig. 13. Effect of change in feed diameter: fifth series. 


and graphs of V,,, against feed diameter con- 
structed. By interpolation from the graphs it was 
possible to calculate accurate values of G V,,, for 
each set of cyclone conditions. Graphs were then 
drawn showing the relationship between G V,,, 
and feed diameter, with all other cyclone vari- 
ables constant for each curve. 

Fig. 14, which is confirmed by many similar 
figures not presented in this report, shows that 
for any given set of cyclone dimensions and feed 
pressure there exists a feed diameter which gives 
a maximum value of GV, and that this feed 
diameter is independent of feed pressure. 

Fig. 15 shows that a very simple linear relation- 
ship exists between this optimum feed diameter 
and the overflow diameter. 

Since G is the volume (or mass) rate of addition 
of water to the cyclone, it follows that G V ,,, is 
a measure of the rate of injection of momentum 
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into the cyclone i.e. momentum added per unit 
’ time. It appears therefore, that for any selected 
; feed pressure and combination of cyclone dimen- 
sions, there is a feed diameter (D,) which results 


5 
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Fig. 14. Average feed velocities against feed diameter and 
calculated G V 4, values against feed diameter. 
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Fig. 15. Relationship between optimum feed and overflow 
diameters. 


in an optimum rate of injection of momentum, 
and that its value bears a simple relationship to 
Within the 
limits of accuracy of the experiments, D, is 


the diameter of the overflow orifice. 


independent of the actual value of the feed 
pressure and also independent of the underflow 
diameter. 

In unpublished data available to the author, 
and relating to dust extraction cyclones, reference 
to the fact that the 
velocity of entry of gases into cyclones can vary 


has been made average 
from a fraction of, to several times, the tangential 


velocity of the gas at the same radius. 

In addition it has been found that maximum 
particle elimination efficiency is obtained when 
the 
Assuming that the same phenomenon exists with 


two velocities are approximately equal. 
hydraulic cyclones and, taking into consideration 
the fact that envelopes of constant tangential 
velocity in the conical portion of a hydraulic 
cyclone are concentric cylinders [3], it is highly 
probable that D, is that feed diameter which 
results in the average velocity at the average 
feed radius (measured from the axis of the cyclone) 
being approximately the same as the tangential 


velocity at the same radius. 


Under these conditions the difference between 
the average velocity of the new feed water and the 
velocity of the water already rotating within the 
cyclone at the same radius would be small and 
the 
minimum. Similarly the rate at which momentum 


resultant turbulent mixing would be at a 


would be “ wasted” in producing undesirable 


turbulence would also be a minimum. 


Such turbulent mixing of the water in the feed 
portion of a cyclone must also result in a similar 
mixing of the solid particles present in the cyclone 


feed. 


The Effect of Change in Shape of the Feed Opening 
Two main reasons offered possible explanations 
of the effect of change in feed diameter on particle 
elimination efficiency (all other variables kept 
constant) : 


(a) The effect of turbulent mixing in the feed 


section of the cyclone as described above. 
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Table 2 


Comparison of volume flows for circular and rectangular feed openings 


| 
Rectangular | Total flows | Underflows % U.F. 
— Circular | Area — 
Length | Width Diameter | Sq. ins. | Rect. Circular Rect. Circular Rect. Circular 
| 4” | 487 15-74 212 | 229 143 | «(145 
| i” | 196 | 13-76 1407 | #183 1-92 13-3 13-7 
590” 12-38 1245 | 154 | 1-56 125 | 125 
B93" 049 10-29 1036 | 116 | 1-18 13 


(b) The fact that all particles did not enter the 
With 


large feed openings some particles must enter 


cyclone at exactly the same radius. 


at appreciably smaller radii than others and 
have a greater distance to settle outwards 
into the which 
reject them to the underflow. 

In order to investigate (b) rectangular feed 


vertical downward flows 


openings were constructed which had the same 
The 


rectangular side parallel to the axis of the cyclone 


area as selected circular feed openings. 
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Fig. 16. Effect of shape of feed inlet on weight discharged 
to underflow. 


was kept as large as possible so that the side 
parallel to any radius of the cyclone was short. 
Duplicate standard tests were carried out using 
solids injected into the feed with circular and 
equivalent area rectangular feed openings. 

In all cases the feed pressures were kept con- 
stant at 15 lbs./sq. in. or 6 lbs./sq.in. Vortex 
finders were 2in. long and overflow diameters 
} in. Underflow diameters were } in. in all cases. 

Both Figs. 16 and 17 indicate that rectangular 
feed openings give slightly improved efficiencies 
for all particle sizes, and that the improvement 


90} 


60 
> 
40 —- 
= CIRCULAR FEED 
© = RECTANGULAR FEED 
Oi>isa 
‘ge T 
| 
| | | 
| | 
° | 
i i } i 5 
FEED DIAMETER 


Fig. 17. Effect of shape of feed inlet on weight discharged 
to underflow. 
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is greatest for the largest feed openings for the 
same feed pressure. Volume flows are slightly 
decreased by changing to rectangular feeds. 


Crude Pressure Drop Studies 


Using a vortex finder 4 in. long, overflow diameter 
} in. and underflow diameter } in., total volume 
flows were measured at a pressure of 15 Ibs. /sq. in. 
for a range of feed diameters. 

In turn the feed jets were unscrewed from the 
cyclone and measurements made of the feed 
pressures necessary to give the same total volume 
flows as those obtained with the various feed jets 
when installed in the cyclone. The pipe line and 
position of the pressure gauge were unchanged. 

Under these conditions the average water 
velocities in the pipe line and through the jets 
were the same with and without the cyclone 
connected. Consequently friction losses were the 
same and the difference between the two pressures 
for each feed jet gave a measure of the pressure 
drop within the cyclone during normal operation. 

In Fig. 18 the results are expressed as the 
fraction of the total feed pressure which was 
“ dropped ” within the cyclone graphed against 
feed diameter. 
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Fig. 18. Fraction of feed pressure dropped in the cyclone 
against feed diameter. 


The results indicate that with large feed open- 
ings a major fraction of the pressure drop occurs 
within the cyclone. This can be explained as 
suggested previously by the excessive turbulence 


which is caused in the feed section of the cyclone. 
For small feed openings the greater fraction of 
the pressure drop occurs in the feed pipe from 
the pressure gauge to the point of entry into the 
cyclone. 

Optimum particle elimination efficiency is 
obtained when approximately 50% of the total 
pressure drop occurs within the cyclone. 


Conclusions 

(1) Using a 3 in. cyclone and a wide range of 
conditions it appears that a }in. feed diameter 
gives optimum solid elimination efficiency for all 
particle sizes. 

(2) To a first approximation optimum solid 
elimination efficiency is obtained under conditions 
which result in the maximum rate of injection of 
momentum for a given feed pressure. 

(3) Long, narrow rectangular feed openings 
give slightly improved efficiencies when compared 
with circular openings of the same cross sectional 
areas, 

(4) It is highly probable that with hydraulic 
cyclones smaller or greater than the 3 in. size used 
in this study, the conditions required to give 
maximum solid elimination efficiencies can be 
obtained by volume flow measurements to deter- 
mine the conditions necessary for maximum rate 
of injection of momentum. 

(5) It is of interest to compare the rough 
economics of inefficient solid elimination using 
a jin. feed with the highly efficient elimination 
obtained using a } in. feed. 

Taking the results given in Table 2 as typical 
examples and remembering that, for the same 
feed pressure and identical cyclone dimensions, 
the jin. feed gives solid elimination efficiencies 
which are considerably greater than those given 
by a jin. feed (Figs. 9, 10, 11, 12 and 13), it will 
be seen that, in spite of the fact that the feed 
area of the } in. opening is only approximately 
one sixth of the % in. feed, the total flow is 10-36 
litres/min. compared with 15-74. 

Consequently three cyclones with } in. feeds 
will treat the same total volume as two cyclones, 
of the same size, with jin. feeds operating at 
the same feed pressure. 
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Pumping costs may be regarded as a direct 
function of the volume handled and will be almost 
identical. For the additional cost of a cheap 
cyclone assembly the same total volume can be 
treated with greatly improved solid elimination 
efficiency. 


The Effect of Change in Overflow Diameter 
Four series of experiments were carried out. Par- 
ticle elimination efficiencies were measured using 
the standard technique. 

In the first series, at a feed pressure of 6 lbs. /sq. 
in., a } in. feed diameter and } in. underflow dia- 
meter were used. Vortex finders were 4 in. long 
and overflow diameters of 4, inches 
were used. 

The second series was carried out under the 
suine conditions except that the feed diameter 
was in. 

For the third and fourth series feed diameters 
were }in. (shown to give optimum efficiencies), 
underflow diameters 4 in. and all vortex finders 
were 2in. long. The conditions for the fourth 
series were the same as for the third except that 


4. 


$! 
— 


89 microns 
80) 
y 


| PARTICLE 
ze 


2 
al ™ 490 
Teli 
10 
20 
22.0 
° 


i i i 
! i 
OVERFLOW DIAMETER 


Fig. 19. Effect of change in overflow diameter : first series. 


the feed pressure was increased from 8 lbs./sq. in. 
to 20 lbs. /sq. in. 

Figs. 19, 20, 21 and 22 summarize the results 
graphically. Corrected weight per cents to under- 
flow were calculated and used for the graphs. 
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Fig. 20. Effect of change in overflow diameter: second 


series. 
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Fig. 21. Effect of change in overflow diameter : third series 
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As the overflow diameter is decreased the 
efficiency of elimination of the finer sizes increases 
steadily but the coarser sizes behave differently 
in that the efficiency increases to a maximum 
and then decreases (appreciably under some 
conditions) with further decrease in overflow 
diameter. 
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Fig. 22. Effect of change in overflow diameter : fourth 
series. 


For a }in. feed, a % in. overflow gives the 
maximum efficiency (all other cyclone variables 
kept constant). For j in. and } in. feed diameters 
the optimum overflow diameter is approximately 
} in. 

Anticipating the results reported in the next 
section it is considered that this fall off in efficiency 
of the elimination of the coarser sizes with further 
decrease in overflow diameter is a direct result of 
the short circuit flow down the outside wall of 
the vortex finder. Referring to a previous paper 
by the author [3] it will be recalled that the 
short circuit flow turns sharply along the bottom 
face of the vortex finder wall to join the main 
overflow stream and results in high inward radial 
velocities. 

In the same paper it was pointed out that, for 
a wide range of conditions, the maximum tan- 
gential velocity was attained at approximately 
0-2 in. from the axis of the cyclone. The position 
of this maximum was independent of the overflow 
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diameter, and from this position to the air/water 
interface the tangential velocity decreased as 
radius decreased according to the relationship 
V ar. 


2 
The centrifugal acceleration ai therefore tends 
r 


to vary directly as r. This relationship holds for 
all horizontal levels beneath the bottom of the 
vortex finder, including the position in very close 
proximity to the bottom face of the vortex 
finder, i.e. in the region where the short circuit 
flow turns beneath the bottom of the vortex 
finder. 

For radii greater than 0-2 in. the relationship 


Vr" held and consequently = increased with 
r 


decrease in radius. 

With large overflow diameters, such that the 
overflow radius plus vortex finder wall thickness 
(0-125 in.), i.e. the radius of the outside wall of 
the vortex finder, is greater than 0-2in., any 
decrease in overflow diameter for constant wall 
thickness results in solid particles in the short 
circuit flow having to pass through regions in 
which the centrifugal acceleration is greater. 
Consequently more particles of a given size are 
eliminated from the short circuit flow. 

However, once the outside wall radius of the 
vortex finder is decreased to less than approxi- 
mately 0-2 in., further decrease results in regions 
of smaller centrifugal accelerations through which 
particles carried by the short circuit flow must 
pass. This results in progressively less efficient 
elimination of particles of a given size carried by 
the short circuit flow, as the overflow diameter 
is decreased. (Constant wall thickness assumed.) 

It has been demonstrated previously [3] that 
the short circuit flow along the outside wall of 
the vortex finder occurs in a region of poor solid 
elimination efficiency. Therefore, only a small 
fraction of the finer particles can be eliminated 
from this flow and, in any case, most of the finer 
sizes discharged through the overflow originate 
from regions in the bottom of the conical portion 
of the cyclone. 

The gain in efficiency of elimination of the main 
fraction of particles in a given fine size range, due 
to the decrease in overflow diameter, will mask 
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the decrease in efficiency of elimination of the veamneas 
same particles from the short circuit flow, size 
Practically all the coarse particles discharged 49:0 microns 
through the overflow are due to the short circuit a 
flow and consequently, as the efficiency of elim- | 80 
ination from the short circuit flow decreases . 
as the overflow diameter is decreased beyond the © 40 . 
optimum size, then the percentage of such xe ys 
particles discharged through the underflow de- 
creases and the overall efficiency for these 
particles also decreases. 
6 20 
The Effect of Change in Vortex Finder Length islilale 
Three series of experiments were carried out in ° l 
which the length of the vortex finder was varied of! 2 3 4 in 
from jin., to 2in. to 4in., measured from the VORTEX FINDER LENGTH comme 
top of the cyclone. Fig. 24. Effect of vortex finder length : second series. 
The conditions were : 
Ist Series : 6 Ibs. /sq. in., jin. feed diameter, } in. | | re 
diameter, underflow 100 49-0 microns 
2nd Series : 15 Ib./sq. in., } in. feed, } in. over- 
flow, 4 in. underflow. 
8rd Series : 15 lb./sq. in., jin. feed, fin. over-| © 60 
flow, in. underflow. 
Inspection of Figs. 23, 24 and 25 shows that an gr | ge 
increase in efficiency is obtained for the finer 2% 
sizes (i.e. the lower settling rate particles) as the 5 20 
length of the vortex finder is decreased. 
° 2 3 4 in 
100 89-0 microns — VORTEX FINDER LENCTH 
>| ° 
5 3 | 10-72 |1-61/ 12-33 
4 [10-98] 1-57} 12°55 
5 40 Fig. 25. Effect of vortex finder length : third series. 
The efficiency of elimination of coarsest sizes 
20 decreases noticeably, up to 10% in some cases, 
with decrease in vortex finder length. 
olL_t Fig. 24 indicates that for some particles the 
i efficiency increases to a maximum and then 
VORTEX FINDER LENGTH =o decreases with further reduction in length of 
Fig. 23. Effect of vortex finder length : first series. vortex finder. 
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The table of flow rates provided with Fig. 25, 
as a typical example, shows that these changes 
in efficiency cannot be attributed to changes in 
throughput or the slight changes in overflow flow 
rate. 

Relatively simple explanations are possible. 

The majority of the particles in the finer size 
ranges which are discharged through the overflow 
have passed down the cyclone outside the zero 
vertical velocity envelope [8] and, in the lower 
regions of the cyclone, have been carried in- 
wards by the water radial velocity components 
to join the high velocity upward flows near the 
central air column. As the radial velocities at 
any radius tend to decrease as the bottom of the 
vortex finder is approached, whereas the centri- 
fugal acceleration at the same radius remains 
unchanged, such particles tend to move outwards 
again. 

With a short vortex finder the time available 
for them to leave the vertical flows, which 
eventually reach the vortex finder direct, is much 
greater than when a long vortex finder is used. (It 
is considered that velocity profiles below the 
bottom of the vortex finder are not appreciably 
changed by the length of the vortex finder.) Con- 
sequently a greater fraction of the particles may 
join the recirculation [8] to the feed section of 
the cyclone, where mixing occurs, and may make 
repeated passes through the main cyclone 
flows. The expected overall result is that shorter 
vortex finders will lead to higher efficiencies for 
the finer sizes. This result is demonstrated in 
Figs. 23, 24 and 25. 

If the assumption that practically all the 
coarser particles discharged through the overflow 
originate from the short circuit flow is correct, 
increase in the length of the vortex finder would 
permit a longer time interval for such particles 
to be slowly eliminated from the short circuit 
flow, and, if the same or similar velocity profiles 
are sustained within the short circuit flow for the 
whole length of the vortex finder, the overall 
result would be a higher efficiency of elimination 
through the underflow, since particles leaving the 
short circuit would be redistributed within the 
feed section and a major fraction would enter the 
main downward flows within the cyclone. 


The Effect of Change in Underflow Diameter 


Since in practical cyclone operation the main 
factor governing the selection of the underflow 
diameter is the quantity of solids to be discharged 
through this orifice, and the total time available 
for the whole investigation had to be limited, only 
two tests were carried out. 

For both tests pressures were maintained at 
6 Ibs./sq. in. In the first a jin. feed was used 
with @ in. overflow diameter and 4 in. long vortex 
finder. Underflow diameters were } in. and } in. 

The second test involved use of a @ in. feed and 
other conditions were the same as for the first 


Table 3 
Flow Rates 
Under-| feed (litres /min.) feed 
flow ——|- -—-- -—- 
diam.| O.F. | U.F. | Total | O.F. | U.F. | Total 
4 12-56 4°36 16-92 9-61 2-97 12-58 


i” | 15-43 | 0-56 | 15-99 | 11-43 


| 11-71 


Fig. 26 summarizes the results. 
Under the conditions of the experiments, i.e. 
constant feed pressure and very low pulp densities, 
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Fig. 26. Particle elimination for different underflow dia- 
meters. 
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it is seen that decrease in underflow diameter 
results in a decrease in efficiency of elimination 
of particles of all sizes. 

Table 4 shows that with small underflows the 
volume flowrate through the overflow increases 
appreciably but that total flows are only slightly 
changed. 

It can be postulated from a general knowledge 
of velocity profiles within the cyclone [3] that 
under the above experimental conditions tangen- 
tial velocity profiles would be very similar but 
radial and vertical velocities in regions beneath 
the bottom of the vortex finder would be appre- 
ciably greater when using smaller underflow 
diameters. Consequently, coarser particles would 
tend to be discharged through the overflow in 
larger quantities. However, much testwork is 
needed to confirm or disprove this crude theory. 


Annulus Experiments 


In order to obtain a quantitative assessment of 
the value of the modification (U.K. Patent 
Application No. 1358 /52), described in a previous 
paper [3], for diminishing the effect of the short 
circuit flow, a simple experiment was under- 
taken. 

The cyclone was set up with a q in. feed, } in. 
underflow and } in. overflow. The vortex finder 
was 2in. long and had a wall jin. thick. An 
annular take off was incorporated around the root 
of the vortex finder, in the flat top of the cyclone, 
with a diameter of }3 in. The construction was 
such that a separate discharge could be taken off 
through the ,), in. wide annulus by a pipe fitted 
with a valve. 

Table 4 


Flow rates for Annulus Test 


litres /min. 


Pressure - 
lbs. /sq. in. Overflow | Annulus | Underflow 
| 
5-9 5-63 3-61 | 0-83 


Feed pressure was kept constant at 5-9 lb. /sq. 
in., and by adjustment of the valve a suitable 
fraction of the total flow was taken through the 


annulus. The usual procedure with batches of sized 
Perspex spheres was adopted but, in this case, 
particles discharged by the overflow, annulus 
and underflow were all collected separately for 
drying and weighing. 

In Fig. 27 the actual weight percents discharged 
through the respective openings are graphed 
against particle size. 


DERFLOW 

3 — 

a 

= 

Vv 

> 

= 

ANNULUS 
‘OVERFLOW 


70 to 90 
SIZES microns 


Fig. 27. Elimination of particles through overflow, under- 
flow and annulus. 


For the coarser sizes the weight per cent dis- 
charged through the annulus tends to a limiting 
value of approximately 15%. As it is virtually 
impossible for very coarse particles to reach the 
annulus from the lower regions of the cyclone, 
it is considered that this phenomenon indicates 
that, of the water discharged through the annulus, 
an amount equal to 15% of the total water flow 
through the cyclone is short circuited directly 
from the feed to the annulus and that, for all 
solids independent of size, a fraction equal to 
15% of the solids in the feed is short circuited 
from the feed to the annulus. 

In order to obtain some information about the 
fraction of all solids discharged through the 
annulus but originating from the lower regions of 
the cyclone, and carried to regions above the 
bottom of the vortex finder by the recirculation 
to which reference has already been made, the 
actual weight per cents discharged through the 
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three openings were corrected in the following 
manner. 

Let 2% of the total volume flow pass through 
the underflow. Then by a previous argument 
(Section 8) (100-7) is the corrected feed and 
(y —z) is the fraction eliminated through the 
underflow by the cyclone action (y% is the experi- 
mentally determined fraction eliminated through 
the underflow). 

A further correction for the feed solids is 
needed for the 15%, short circuit i.e. (100 — z—15) 
is the final feed corrected for both underflow and 
short circuit to the annulus. 

If z is the actual weight %, of selected solids 
discharged through the annulus then (z — 15) is 
the corrected fraction. 

On this basis, with double correction, the 
corrected weight per cent to the underflow is 


(y — 2) 


(100 —z — 15) x 100, the corrected weight per 


cent to the annulus — 
(100 


the remainder to the overflow. 

Fig. 28 is of considerable interest. Remembering 
that the short circuit direct from the feed has 
been eliminated mathematically by the above 
corrections, inspection shows that the coarse 


= = 


WT % TO RESPECTIVE DISCHARGE (CORRECTED FOR UF AND ANNULUS) 


° 10 so 80 70 


microns 


Fig. 28. Particle elimination through overflow, underflow 
and annulus corrected for annulus and underflow factors. 


particles have a very high probability of discharge 
through the underflow as a result of cyclone 
action, but very low probabilities of discharge 
through the overflow or the annulus. 

Very fine particles have a high probability 
of discharge through the overflow but low 
probabilities of discharge through the underflow 
or annulus. 

Under the conditions of the experiment part- 
icles of intermediate size (i.e. in the region of 


OVERFLOW 
BN 
\ 
\P 
OVERFLOW 
DIAMETER 
UNDERFLOW 


Fig. 29a. Vertical section through cyclone showing con- 


struction 


30 microns) have low probabilities of discharge 
through either the overflow or the underflow. They 
have, however, a relatively high probability of 
discharge through the annulus. 

Flow patterns within the cyclone must be 
affected by withdrawal of a portion of the flow 
through the annulus, but it is considered highly 
probable that the result presented in Fig. 28 
provides a qualitative picture of the behaviour of 
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intermediate particles during normal cyclone 
operation. As confirmation it may be stated that 
during the extensive testwork described in this 
report it was often noted that some particles 
required excessively long periods of time before 
they were completely eliminated from the cyclone, 
and the recirculation, already described, was 
easily visible in regions above the bottom of the 
vortex finder. 


CONSTANT HEAD 
FEED CONE 
GAUGE 
CYCLONE 
j 
TO WASTE 


UNDERFLOW OVERFLOW 


CONSTANT HEAD 
FEED CONE 


j 


TO WASTE UNDERFLOW OVERFLOW 


Diagrammatic arrangement of cyclone circuit 
for high and low pressure operations. 


Fig. 29b. 


The main significance of this experiment is that, 
for a given set of cyclone conditions, some parti- 
cles have a small probability of discharge due to 
the cyclone operation (i.e. centrifugal forces 
within the cyclone) through both underflow and 
overflow. Such particles tend to recirculate 
within the cyclone and are gradually eliminated 
by two methods. After remixing in the feed 
section of the cyclone a fraction joins the down- 
ward flow in which, as described in Section 8, all 


particles independent of size have a probability 
of approximately one of being discharged through 
the underflow. A further fraction joins the short 
circuit flow down the outside wall of the vortex 
finder and a portion is eliminated through the 
overflow. Some particles will be thrown out of 
the short circuit flow to rejoin the circulation. 

This phenomenon must be of considerable 
importance during continuous cyclone operation 
since particles of a certain size range will tend to 
accumulate within the cyclone and result in an 
increase of pulp density and effective viscosity. 
Methods of overcoming this effect are under 
consideration. 


5. CONCLUSIONS 
For the conditions, under which the tests reported 


were carried out, (i.e. almost infinite dilution) the 
following statements can be made : 


(1) Fora given set of cyclone dimensions there 
exists a power relationship between d,, and feed 
pressure (or throughput) of the form 

1 1 

(where m and n < 1). 

The available evidence indicates that the values 
of m and n are dependent on cyclone dimensions. 


(2) For optimum solid elimination efficiency 
with a 8 in. cyclone, the feed diameter should be 
jin. To a first approximation this feed gives 
optimum efficiency for all cyclone conditions, and 
corresponds to the conditions under which 
maximum rate of injection of momentum is 


obtained. 


(3) Long, narrow rectangular feed openings 
give slightly improved efficiencies when com- 
pared with circular openings of the same cross- 
sectional areas. 

(4) With large feed openings (e.g. ¢ in.) 
excessive turbulence in the feed section of the 
cyclone results in large pressure drops within the 
cyclone. 

(5) Reduction in overflow diameter, with all 
other cyclone variables constant, results in 
increased efficiency of elimination of the finer 
sizes, but for coarser sizes the efficiency increases 


271 


. 
L . 
953 
GE 
Pos 
CYCLONE 
VALVE 
PUMP 
| 
= 
| 


DD. F. Keusaus 
to a maximum and decreases with further 
decrease in overflow diameter. This phenomenon 
is considered to be a direct result of the short 
circuit flow down the outside wall of the vortex 
finder. 


(6) As the length of the vortex finder is 
decreased an increase in elimination efficiency is 
obtained for the finer sizes but appreciable 
decreases in efficiency are obtained for the 
coarser sizes. The latter may be attributed to 
the fact that practically all the coarser particles 
discharged through the overflow originate from 
the short circuit flow. 


(7) Decrease in underflow diameter for the 
same feed pressure, with all other variables 
maintained constant, results in a decrease in 
efficiency of elimination of all particles in spite 
of the fact that only small decreases occur in 
total volume throughput. 


(8) Under conditions of normal cyclone opera- 
tion, particles in intermediate size ranges tend to 
have low probabilities of discharge through both 
underflow and overflow as a result of cyclone 
action. Such particles tend to recirculate from 
the lower regions of the cyclone to rejoin the 
mixing section near the feed opening. They leave 
the cyclone either by joining the short circuit 
flow and then leaving via the overflow, or by 
being carried down by the strong vertical flows 
near the conical cyclone wall in which all particles, 
independent of size, have a high probability of 
leaving via the underflow. 


(9) From (1), (2), (5) and (6) above it is 
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considered that no simple power relationships 
can be derived to describe the change in solid 
elimination efficiency resulting from alteration 
of feed opening, throughput or overflow diameter, 
to cover the normal operating range which must 
include the most efficient conditions. 

Two main phenomena are responsible for the 
complexity of the relationships, namely, the 
variable degree of turbulent mixing due to 
“shock” effects in the feed portion of the 
cyclone, and the importance of the short circuit 
flow down the outside wall of the vortex finder. 


(10) Finally, from (5) and (6) above, neither 
the d,, method for assessment of efficiency nor 
the method used by Fircn and Jounson can be 
used to describe the overall efficiency for a com- 
plete size range with any degree of significance. 
The only reliable method involves a comparison 
of efficiency curves and due consideration of the 
metallurgical requirements demanded of the 


separation. 
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Summary— Dimensional analysis in the case of fluidization, leads to expressions which apply 
equally well to fixed beds, to fluidized beds or to the “ critical conditions in the bed, when 
fluidization starts. 

For porous beds, the ratio between the pressure drop and the mass velocity of the fluid expresses 
the resistance of the bed. This leads to the definition of the * conductivity ” of the porous 
medium, which depends on the physical properties of the fluid, the shape of the particles and 


their geometrical arrangement. The latter can be represented by a function of « the porosity 


(or voidage fraction) which is called the “ vacuity ” g : 


For a given fluid and set of particles, the conductivity can always be represented by a function 


of the vacuity, Kgs for a fixed bed, and Kes", for a fluidized bed. 
Comparison between Kzy and Kgs" leads to an expression for v, the average velocity of trans- 


lation of the solid particles in the fluidized bed. 
Very few measurements are required to gather data for the control of conditions in porous 


beds, either fixed or fluidized. 


Résumé—Appliquée a la fluidisation, analyse dimensionnelle permet d’établir des relations 
qualitatives et quantitatives qui s‘appliquent aussi bien & toutes les couches poreuses, que les 
particules soient immobiles ou soient fluidisées et également dans cet état critique particulier 


ou s‘amorce la fluidisation. 
L’auteur montre que, pour toutes les couches poreuses, le quotient de la perte de charge 
par la densité de courant fluide représente la résistance de la couche & 'écoulement. Ceci con- 
duit a définir une “conductivité” d'un milieu poreux qui dépend a la fois des propriétés du 
fluide, de la forme des grains et de leur répartition géométrique. 
Cette dernitre peut étre représentée par une fonction particulitre du “ degré de vide 


ou porosité « que l'auteur propose d’appeler vacuilé g : 


Pour un fluide et un ensemble de grains donnés, la conductivité peut toujours s'’exprimer 
par une fonction simple de la vacuité, Kzgy pour une couche immobile et Kz" pour une couche 
fluidisée. 

De la comparaison entre Kzs et Kzj", on aboutit a une expression donnant v, vitesse moyenne 
de translation dans le lit fluidisé. 

L’auteur montre comment, a l'aide de quelques déterminations expérimentales sur un 
échantillon de particules solides, on peut obtenir toutes les indications numériques suffisantes 
pour permettre d'estimer les conditions actucllement réalisées dans écoulement au travers de 


la couche poreuse, qu'elle soit immobile ou fluidisée. 


Amongst the numerous unit operations which are | fluidized beds, we deal with an intimate mixture 
the proper tools of the chemical engineers to meet | of solid particles and fluid flowing at such a velo- 
the requirements of the chemical industry, I think | city that the particles are lifted and supported by 
that fluidization deserves special attention. In | the fluid. In the particular state of fluidization, 


* This was delivered as a lecture to the University of London, at University College in January, 1953. 
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the particles can be considered as free moving 
bodies : inside the bed, each of them follows an 
irregular course according to the chances of 
momentum exchanges with other particles or 
with the fluid. Through these repeated momentum 


exchanges, we achieve a very important reduction 


of the normal resistance to the heat and mass 
transfer between the solid and the fluid. 
Fluidization is now one of the most useful 


methods of achieving exact control of all the 
factors regulating any chemical exchange between 
fluids and solids. The details of any of its many 
industrial applications remain outside the scope 
of this paper. But I consider that the development 
of the theory of fluidization gives a very striking 
example of the methods special to chemical 
engineering science for solving its own problems, 
Today, we know far more about the essentials of 
the theory of fluidization than we know about some 
very old unit operations, grinding for instance. 
Like any scientist, before any attempt towards 
theory, the chemical engineer must start some 
experiments. 
relation between the flow of some particular fluid 
bed of solid 
He can measure u the linear velocity 


and the pressure drop through a 
particles. 
of the fluid or preferably the mass velocity, G, 
mass of fluid crossing in a unit of time the unit 


of section of the empty column. Usually called 


“mass velocity,” G is a current density and its | 


relation to the linear velocity is given by 


G=up G=Me"L? 


The laws governing the flow of fluids through 
porous beds are well known and between very 
large limits the pressure drop AP remains pro- 
When G 


becomes very high, the pressure drop increases 


portional to the mass velocity (Fig. 1). 


more slowly than the mass velocity. Finally, the | 


mass velocity can continue to increase at a 


practically constant pressure drop: we have a 
fluidized bed. But if we stop the flow and start 
again our experiments, we observe new values 
for G or AP: they vary according to the con- 
ditions of our bed. 

To give a quantitative representation of these 
conditions of our porous bed, we will call ** packing 


fraction +r (in French, tassement the ratio 


He will determine for instance the 


Pressure drop and rate of flow (spheres of uniform 


Fig. 1. 
diameter) 


of the occupied volume to the apparent volume, 
We will call * voidage fraction ” or porosity « the 


We 


are also interested in the specific height of our 


ratio of the voids to the apparent volume. 


bed which can be expressed as the height of a 
unit of mass of particles in a column of unit 


area. The specific height A is oby iously related 


to the apparent volume and we have : 7A 1 /p’; 


e+fT a 


Dealing for instance with regular spheres of 
constant diameter, we know that the porosity will 
vary according to pure geometrical laws, from a 


minimum value 0-259 for a rhombohedral distri- 


bution to a maximum 0-476 for a cubic distribu- 
| tion. It is the 
specific height of our bed at its maximum porosity 


convenient to take as a unit 


To obtain reproducible determinations in our 
fluidization experiments, we must start with a 


bed of maximum porosity, « For beds with 


closer packing, before fluidization starts, we have 


some local rearrangements of our spheres, and a 


corresponding increase of the free area for the 
flow with a small decrease in the pressure drop. 
With proper care, we observe a well-pronounced 
break in our curves : fluidization starts just above 
a threshold and requires a critical value G, for the 
mass velocity. 

Before this the 
pressure drop increases, as the increased velocity 


threshold, when G grows, 


increases the friction against the walls of the 


small irregular channels of our bed. Inversely, 
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through the friction, every sphere suffers an 
increased vertical pressure, opposed to its weight. 
For a certain critical velocity, G,, the friction 
forces exactly balance the gravity : our spheres 
remain still tangent together, but they cease to 
support each other. For a very small increase 
above G,, the spheres are separated and we are 
in the fluidized condition. 

We know that the friction forces will be a 
certain function of the diameter of our spheres, 
will be proportional to the linear velocity, and 
will depend on the viscosity of our fluid. These 
friction forces will be equal to the weight of our 
particles, which is proportional to some function 
of their diameter, their buoyancy and g. We can 
write 

xd (p' — g = uy 


x am (p’ p)” pg” 


uy p G 


From the dimensional analysis of the critical 


conditions before fluidization starts, we can 


obtain the values for the exponents m, n and z : 
m 2 n 


G,, will be proportional to the buoyancy of our 
particles, to the specific gravity of our fluid, to 
the reciprocal of its viscosity and to the square 
of the diameter of our spheres. The height of our 
In fact, at 
the critical point, the same critical velocity will 


bed does not appear in our formula. 


totally support beds of very different heights. But 
the pressure drop will increase proportionally to 
the height. 

Now we will gather the answer from the experi- 
ments. Without going into the details of the 
many determinations published for critical mass 
will give briefly significant 


velocity, we some 


results. Leva {1} and colleagues’ results can be 
summarized by the relation : 

0-005 -—(p 

(1 €q) 


Miiier and Locwinuxk [2] propose the following 
expression 

G, 0.00125 d* (p' p' 


Assuming perfect spheres of constant radius in 


27 


imposed by 


MILLER’s experiments, with 0-476 ; l—e 


0.2065, their relation becomes : 


3 


G, 0.00605 dB (p’ 


For Van HEERDEN [3] and colleagues the critical 
mass velocity is given by : 


3 
G, = 0.0065 


The concordance between the three sets of 
experiments is impressive and we come to the 
conclusion that our dimensional analysis of the 
conditions starting 


existing at the point of 


fluidization is correct. G, varies with the square 


0 
of d, with the expected dependence on the specific 
gravities and viscosity. But we have not yet 
obtained any information on the pressure drop, 
necessary to maintain the mass velocity G,. 
To go further, we will remember that, in the 
critical stage we consider that our particles are 
still in material contact together. Except for the 
amount of pressure between the particles, at the 
critical point, no differences exist in the con- 
ditions of the porous bed and it is legitimate to 
apply the well known laws on the flow of fluids 
Since the 
Kozexy [4] and Carman 


researches of 
that 


PoIseEvILLE’s law applies with some modifications 


through porous bed. 
[5], we know 
the irregular shape of our small 


channels : 


dP 1 
32 
db d* 
u. € u 
a N 
Nor d* a D* «4 
d be 
l 
dP (1 «)* 
pua 
all, 
«= (= ) = 2-467 = 25 
a=-rd? 
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AP (i-? G ¢ 
= 180 
(a) L d* 
AP 
(b) — p)g(1 = 
180 (1 — 
p @ 
(ec) Gy = 0.0055 | p) 


Not only did we find the same form of depend- 
ence between G, and the porosity supported by 
the experiments of Leva and Van Herrven, but 
the numerical value of the constant lies between 
the extreme values found in the experimental 
work. 

To obtain in expression (c) the value of G,, the 
critical mass velocity when fluidization starts, we 
make in (b) the pressure gradient AP/L as given 
in (a) equal to the weight of all the particles 
present in a volume unity : our theory is based 
on the assumption that all these particles are 
supported now by fluid friction. 

But expression (a) applies to all fixed beds of 
porous particles and shows the interdependence 
between mass velocity G, pressure gradient APL 
and the variable porosity « of any bed of definite 
particles. The four different straight lines in 
Fig. 1 correspond to the values taken by « when 
spheres of constant diameter are assembled 


together in regular geometrical arrangement : 


= 0-259 rhombohedral packing. 
€, = 0-802 tetragonal packing. 
0-395 orthorhombic packing. 
0-476 cubic packing. 


To each value of the porosity corresponds a 
value of the specific height of the bed of particles 
(Fig. 2). 


v M 
Packing fraction + 


Porosity « = = 1 
v 


Specific height A 1/p’ r. 


It is worthwhile to transform expression (a) : 


2 a 
P p 


é 30 040 ows 060 | Ors 
c Fluidized 
Fined bed bed 


— — | 


Fig. 2. Specific height of porous beds 
(spheres of uniform diameter). 


Here as well as in all its transfer problems, 
chemical engineering science finds a generalized 
form of Oum’s law : 


Driving force 


Rate of flow Resistance 
Ss 


The quotient of the current density by the 
gradient of the driving force is the reciprocal of a 
resistivity : it is a conductivity. The conductivity 
corresponding to this transfer problem has the 
dimension of time. 

In relation (c), we will note that, at the critical 
point, the current density is identical with the 
product of the bed’s conductivity by the apparent 
weight of all the particles enclosed in a volume 
unity : 


3 
G,=0.0055 . . 


—p)8 (1 


In the mathematical expression of the con- 
ductivity of a porous bed, it is obvious that we 
shall find a term related to the physical properties 
of the moving fluid : this property is the recipro- 
cal of its kinematic viscosity. Leaving apart the 
“ kinematic fluidity ” of the fluid, our expression 
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of the conductivity gives the * permeability ” of 
the bed, a characteristic property depending upon 
the diameter and shape of the particles (terms 
@® ¢*) and upon their geometrical distributions. 
We propose to call * vacuity ” of the bed and to 
represent by the symbol U the term (1 — «)? 
which describes the arrangement, and we can 
write for any porous bed ; 


(kinematic fluidity) (area of 
equiv. sphere) (vacuity). 


conductivity 


GL 
ad- @*)0 
AP (pp) 0-055 =a) 


Fora given fluid, it is clear that the conductivity 
of all porous beds made with the same particles, is 


a linear function of the vacuity of the bed ; 


GL 
AP 


AKO 


From direct measurement of the height L of a 
porous bed of constant area S, and containing a 
mass M of particles in a volume V, we can obtain 
the specific conductivity of the porous bed of 
unit area and specific height A, corresponding 
to a mass unity. The pressure drop takes a new 
numerical value but the general expression still 
remains : 

GA 
AP 


KO 


When the porosity varies between 0 and 1, the 
take 
infinity. Table 1 gives the corresponding values 
of «, U, and log U. 

In any porous bed of constant area S, it is 


vacuity can any value between 0 and 


always possible to measure L and G in various 
conditions of the porous bed : fixed bed or fluid- 
ized bed. From L, it is easy to obtain A, then 
« = 1 —1/p’ A and finally the present value of 
U. Table 1 gives a series of values for K A and 
log K A, K being p’ (1 
to 0-476, its normal value for a bed of spheres of 


- €&) and «, taken equal 


constant diameter. 

Dealing with fluidized beds, when AP remains 
constant and equal to AP, (the value at the 
critical point), Leva has been able to show that 
the conductivity grows with the nth power of 
the vacuity. The experimental values of n vary 


between 2-5 and 3-5, according to the size and 
shape of particles ; ” remains independent of the 
nature of the fluid. 
G, A 

f K 
AP, 
‘expanded bed as 
In the 


conditions of our expanded bed, the specific 


Now we can imagine an 
opposed to a fluidized bed. idealized ” 
height will be the same, the average value of the 
mean distance between our particles will remain 
the same as in the fluidized state but our particles 
will remain immobile : exactly like in a porous 
bed, the conditions of an expanded bed can be 
represented by a straight line of slope unity. 

For every value of the vacuity, we have a 
definite value for G, and for G,, the mass velocities 
corresponding to this fluidized and “ idealized ” 
conditions. Leva calls fluidization efliciency the 


value of the ratio. 


This 


intensity of the agitation in fluidized beds. 


ratio is very useful to appreciate the 
We will remember that Pressure multiplied by 
Volume gives Work. The product of the pressure 
drop by the total rate of flow through the entire 
bed gives the power required to maintain the 
bed in the fluidized or in the “ expanded ” state. 


Fluidized conditions : 


= AP. D*/4. G,/p 


Expanded state : 
AP. D*®/4 . G,/p 
obtain 


From the difference, we the 


required to maintain the agitation of our particles 


power 


in the fluidized bed. 


W, = W,—W,= AP. G,. 


Such a power is consumed by friction, each 
individual with 
velocity v hitting many times the other bodies. 
Calling F a friction factor we have : 


particle moving an average 


Wo = Fgv lm 


a 
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where, for the critical conditions, we have : 


AP 

om =p 4 : (p P) 
Ag 

and remembering that for fluidized bed AP 

constant = AP... we obtain the value of the mean 

velocity of our particles in a fluidized bed. 


v=G,.7 

In this way, we have developed some kind of 
kinetic theory of our fluidized state and we are 
now in a position to attempt some analysis of the 
conditions of heat transfer in a fluidized bed. 

It is classical in heat transfer problems to refer 
to A, the heat transfer coctlicient, the dimensions 
of which are M@*T': h depends from both 
the properties of the system and the conditions 
of flow. 

The properties of the system governing the 
heat transfer will probably be its thermal con- 
ductiv ity k, its viscosity ja, its specific mass p and 
its heat capacity C, and the conditions of flow 
which shall be represented by the value of v, the 
average velocity of our particles, and by some 
linear dimensions D. We can now write a general 
expression of h : 


h=« vke Ce 
and the dimensional analysis leads immediately 
to an expression of classical form : 
hD 
x 
k k 


x DG,» 
The experimental determinations of A for 
fluidized beds are not very numerous; their 
results are often contradictory. From his experi- 
ments, Leva considers that g and r can be neg- 
lected and gives for the heat transfer coefficient : 


x 
h kG,» 


We think that the experiments of Leva are 
not suflicient to prove that linear dimensions are 
without influence on the heat transfer coefficient. 
The fluidization efficiency depends for a particular 


bed on the value of the slope of the fluidization 
line and this slope varies with the diameter of 
the particles. Some authors have found explicit 
influence of this diameter but from our dimen- 
sional analysis, we know now that to detect the 
role of linear dimensions, it cannot be sullicient 
to experiment with particles of different diameter 
or shape in the same column, but it is essential 
to vary also the size of the bed. 

It is very convenient to represent on a graph 
(Fig. 3), the relationships between all the charac- 
teristics properties of the flow through our beds. 


on 
| 


Phomsoredra/ 


» bp 
4 
© s >SSs S 
Fig. 3. 


Curve 1 represents the relation between ¢ (ordinate 
axis) and log U (abscissa axis). From curve 1, it 
is useful to transfer the values of the porosity on 
a horizontal line. Curve 2 gives the values of 
log A as a function of log U. 

We have given in Table 1, the series of numerical 
data needed for the construction of curves 1 and 
2. According to our theory, curves 1 and 2, with 
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Table 1. 
Observations 
P KA log KA oO log @ (spheres of uniform diam.) 
Ol 0-585 0-233 0-00123 2-010 
0-2 0-655 0-184 0-0125 1-903 
0-25 0-028 1-156 
0-26 0-70 0-155 0-032 1-495 Rhombohedral, minimum porosity 
0-3 0-748 0-126 0-055 1-260 
0-302 O-75 0-125 0-056 1-252 Tetragonal 
0-35 0-906 
0-395 0-87 0-06 0-166 0-78 Orthorhombic 
0-476 1 0 0-393 0-406 Cubic, maximum porosity 
0-500 1-05 0-021 O-S 0-301 | 
0-55 0-821 0-086 
0-565 120 0-079 1 0 
0-6 1-3 0-114 1-35 0-130 
0-65 2-24 0-35 
0-666 1-57 0-196 2-66 0-425 
O-7 1-75 0-243 3-81 0-581 
0-75 2:1 0-322 6-75 0-829 
0-80 2-6 0-415 12-8 1-107 
O-85 3:49 0-543 1-436 
0-9 5-25 0-720 72-9 1-863 
0-95 10-5 102 343 2°535 


their corresponding scales apply to any kind of 
fluid flow through any kind of porous bed. 

The control of any fluidized process requires 
the knowledge of the conductivity of our bed and 
the average velocity of our particles. Few 
experimental determinations are needed but they 
can be achieved on a small scale, at normal 
temperature, and with air for which accurate 
kinetic fluidity values are readily available. 

For any kind of particle, of variable size and 
shape, it is always possible to measure the 
porosity of the bed when we know the specific 
mass p’. In a column of constant section S, we 
will dispose a mass M of particles and measure 
the height of the bed L : we obtain immediately 
A= LS/M and the packing fraction or the 
porosity: r= 1—e=p' A. 

It is worthwhile to measure various values of 
the porosity for various conditions of packing. 
Dealing with air, we can measure, for every 
conditions of packing, the values of the mass 
velocity G and pressure drop AP. In this way, 
we can obtain a reliable experimental value of 
¢*. For irregular particles, it scems that the 


correct value of the dimaeter of the equivalent 
sphere is more accurately obtained when operating 
the bed at its minimum porosity. 

To obtain the maximum porosity, we will 
increase the rate of flow until fluidization con- 
ditions are well established. By reducing gently 


the rate of flow, we can measure A, and ¢,. As 


0 
explained before, it is convenient to take A), as 
the unity of specific height for our beds. Accord- 
ingly, in Fig. 2, the point corresponding on curve 
2 to log U, (and to «, on the horizontal scale) will 
be taken as the origin of the vertical scale of 
log A. 

When we carry out the determinations of A, 
and «,, we can measure also together G, and 
AP,. It is also convenient to take the same 


point as the origin of the vertical scale for the 


function log This is equivalent to take for 


AP 
respective units G, the critical mass velocity, 
and AP, the pressure drop, together with A,. 
Dealing with mass velocities above G,, can meas- 
ure different values of G 3G,..4G,...6G, 
and obtain the corresponding values of L, A, «, 
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U. In fluidized conditions the pressure drop | we can utilise our graph by converting the critical 


remains constant AP,. It is easy to obtain the 
value of the exponent n in the expression : 
G, L 
Al’, 


K 0" 


equivalent to : 
log G, A = n log + constant 

In Fig. 3, this corresponds to a straight line, 
of slope n, which we start from our chosen 
origin. A straight line of slope 1, corresponds 
below curve 2 to the flow conditions in fixed 
beds and above curve 2, to the “ idealized” 
expanded bed. In Fig. 3, we have taken for «, 
the value 0.476, maximum porosity of a bed of 
spherical particles with constant diameter. 

Once constructed for a definite set of particles, 
for the control of any fluidized bed, the graphs 
of Fig. 3 give immediately all the necessary data 
by the single determination of the pressure drop 
through a known height of the bed. From this, 
we can deduce the present value of the specific 
height. In our graph, curve 2 gives log A 
segment PH; we obtain directly the logarithm 


of the present mass velocity G, (segment HF); 


the present value of log.G, is given by segment 
H F. We obtain easily the fluidization efliciency 


and the average velocity of our particles. 


If instead of working our bed with air at | 


normal temperature, we deal with any other gas 
in different conditions of temperature and pressure 


mass velocity with air at normal temperature to 
the new value of the critical velocity by intro- 
ducing the ratio of the values of the kinematic 
fluidity in the two sets of conditions (we can 
neglect the effect on the buoyancy of the particles 
of the change in specific gravity of the gases). 


I do not pretend that the above theory is 
either complete or definitive : many important 
details have been left out. In our kinetic theory 
of the fluidized bed, we have been compelled to 
accept coarse simplifications. In practice, our 
particles are neither spherical nor uniform in 
size. Many experiments and plenty of work 
remain to be done before achieving a complete 
analysis of the mechanism of fluidization: we 
need for instance better methods of expressing 
the size distribution of our particles ; until now, 
the influence of their shape has been compl tely 


neglected. Nevertheless, we are now able to give 


| with some accuracy the basic data for the design 


| of a fluidized bed converter. 


When we remember that the large utilisation 
of fluidization is not yet ten years old, it is clear 
that theory has not been slow in super. eding 
empiricism. I think that we can recognize a true 


science when we find any organized knowledge 


| helping us to solve some problems by its own 


methods. Fluidization is for me a sound example 
of the special methods which belong to chemical 


enginecring science. 


NOMENCLATURE 


Symbol Dimension Signification 
a L" Specific area of a particle 
Cc Heat capacity of the fluid 
d L Diameter of the equivalent sphere to a particle 
d, L Diameter of channel for flow in porous bed 
D L Diameter of the fixed or fluidized bed 
F 0 Friction factor 
Le? Acceleration of gravity 
G ML“e- Mass velocity of the fluid (empty column) 
h Heat transfer coefficient 
ML@* T-! Thermal conductivity 
L L Present height of the bed 
N Oo Number of channels in the section of the bed 
AP ML-“@° Pressure drop through the bed 
“ Le Linear velocity (empty column) 
v Le Average linear velocity of a particle in fluidized bed 
w ML*@* Power required for maintaining flow conditions 
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Symbol Dimension Signification 
m,n, p,q 0 As exponents 
z,y,% 0 As exponents 
° 0 As exponent for specific proportion of the particle 
0 0 As subscript, characteritic of critical conditions 
of 0 As subscript, characteristic of fluidization 
e 0 As subscript, characteristic of expansion 
a 0 As subscript, characteristic of agitation of bed 
A M"*L® Specific height of the bed in present conditions 
ML“ Fluid viscosity 
Specific gravity 
a Numerical factor 
« 0 Porosity of the bed 
” 0 Fluidization efficiency 
T 0 Packing fraction of the bed 
¢ 0 Shape factor of the particle 
o 0 “ Vacuity of the bed (gy = — 
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Extraction of Acetic acid from water 


2. Ethyl acetate—acetic acid water 


F. H. Garner and S. R. M. Ex.is 


The Chemical Engineering Dept., The University, Birmingham, 15 


( Received 11 November 1953) 


Summary — Solubility and liquid-liquid equilibrium data are given for the system ethyl acetate- 
acetic acid-water at 30°C, 40°C, and 50°C ; the solubility curve and distribution ratio are affected 
only slightly by change in temperature. The predicted distribution coefficients are lower than 
experimental coefficients. Papers are to follow giving the vapour-liquid equilibrium data for 


the system ethyl acetate-water-acetic acid. 


Résumé Présentation pour le systéme acétate d'éthyle acide acétique eau 30°C 40°C 
50°C, des solubilités et des compositions d‘équilibre liquide-liquide ; les variations de tempéra- 
ture agissent trés peu sur les couches de solubilité et sur les coefficients de distribution. Les 
valeurs caleulées pour ces dérniérs sont plus faibles que les valeurs expérimentales. 
Pour le méme systéme, d'autres mémoires suivront pour donner les valeurs relatives a l'équilibre 


liquide — vapeur. 


INTRODUCTION 

EaGLESFIELD, and Snort [14] have 
suggested that for less than 15.0 per cent by 
weight of acetic acid in the feed, ethyl acetate 
is one of the best solvents available. Above 15-0 
per cent by weight of acetic acid in the feed, the 
ethyl acetate may be mixed with a proportion 
of benzene. 

Although the system ethyl acetate-acetic acid- 
water is important, industrially the published 
data is of a limited nature. The mutual solubility 
curve was determined by Mion [15] at 30°C. Tie 
line or liquid-liquid equilibrium results were 
reported by Sononr and WarHADPANDE [16] for 
30°C while and Suort [14] 
give distribution coefficient data for 15-0°C. Roy 
[17] using the apparatus described by Hunter 
and Nasu [5] determined equilibrium data at 
30°C and 40°C. Equilibrium results were also 
obtained by Roy [17] at 50°C and 60°C by using 
a stoppered flask. 

It was thought desirable that the results of 
Roy [17] should be repeated using a modified 
form of the equilibrium apparatus of Smrrx and 
Bonner [18]. 


Purity or Mareriats Usep 
Analytical grade chemicals and a laboratory 


2R2 


source of distilled water were used. Ethyl acetate 
and benzene were fractionated in a thirty plate 
column at a reflux ratio of 25:1. The solvents 
and the acid were tested for purity by refractive 
index and boiling point determinations, 


Substance Refractive Index Boiling Point 
at 20°C 
Water 1.8320 100.0 
Acetic acid 1.3680 118-5 
Ethyl acetate 1.8700 77-1 
Benzene 1-5011 80-1 


EXPERIMENTAL 


The method of analysis of acetic acid and the 
determination of solubility data was as discussed 
in Part 1 benzene-acetic acid-water [19]. 

The tie line data was determined in a modified 
form of the apparatus of Smiru and Bonner [18]. 

Solubility data for ethyl acetate-acetic acid- 
water at 30°, 40°, 50° and 60°C is given on Table 
4. Table 5 gives further solubility results at 30°C 
for the system 85/15, ethyl acetate /benzene- 
acetic acid-water. 

Equilibrium or tie line data for the system 
ethyl acetate-acetic acid-water at 30°, 40° and 
50°C is given in Table 6. 
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F. H. Garner and S. R. M, Exv.is : Extraction of acetic acid from water — 2 
Acetic Acid. 


FIG. 3 


Ethyl Acetate. 


Discussion oF REsuLts 
A comparison of the results given on 
Table 4 shows that the spread of the 
solubility curve decreases only slightly 
with temperature. The equilibrium data 
is almost independent of temperature 
changes. 

Fig. 3 shows the equilibrium data 
at 30°C plotted on triangular graph 
paper. These results and those at 40°C 
and 50°C are satisfactorily correlated by 
Bacuman [6] and Orumer-Tosias [7] 
plots. Fig. 4 plots the equilibrium data 
as a distribution curve and compares the 
results of several investigators. The 
experimental results given here are in 
close agreement with those of Sonont and 
WarRHADPANDE [16] and of HaucuTron 
[20] determined at 25°C. The data of 
EAGLESFIELD, and [14] 
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Table 4 give a higher distribution coefficient in favour 
Ethyl Acetate- Acetic Acid-Water. of ethyl acetate at high concentrations of acctic 
acid, 
For the system 85,15 ethyl acetate, benzene- 
acetic acid-water only a few equilibrium deter- 
| Acetic Acid | Water | Ethyl Acetate | minations have been made. With a 19-5 per cent 
: ec gm gm gm concentration of acetic acid in the solvent layer 
ee there is a concentration of 24-0 per cent acetic 
30 49 | 6-4 88:7 
acid in the aqueous layer. 
12-4 14-0 73-6 
14:8 69-7 
20-8 26-2 53-0 
22-0 43-6 
21:8 48-9 29-3 
58 85-9 8-3 
9-3 80-3 10-4 
21-8 49-6 28-6 
21-7 31-1 47-2 
19-3 22-4 58-3 / 
49 6-5 87-6 
9-4 10-6 80-0 
13-2 14-0 72-8 = 
17-8 210 62-2 
20-2 28-3 
21-2 370 418 
4:3 77 
10-0 80-0 10-0 
13-5 74:2 12:3 
16-5 66-5 
20-2 54-0 25-8 
o-4 10-7 79-9 
| 13-2 14-0 72-8 FIG.S 
144 16-5 69-1 
17-7 21-6 60-6 Table 5 
19-6 27-8 52-6 85/15 Ethyl Acetate / Benzene- Acetic Acid-Water. Solubility 
20-6 37-2 42-2 Data 
5-8 86-3 
9-3 80-9 o-8 Acetic Acid | Water Ethyl Acetate 
141 72-6 13-3 rc gm yum gm 
17-0 65-0 18-0 
19:3 55-0 25-7 9-6 6-6 
. — -- 20-3 | 64-6 
2 60 5-0 6-7 88-3 20-6 30-6 39-8 
150 15-0 70-0 30-6 33-3 35-9 
18-5 26-0 53-5 32-7 38-3 29-0 
18-7 28-0 53-3 36-5 15-0 
48 86-7 8-5 36-0 50-0 140 
7-2 83-7 9-1 29-0 65-7 3-3 
130 739 13-1 20-7 76-5 2-8 
16-3 66-3 17-4 20-0 77° 2-8 
19-0 5655 24:5 0-0 98-3 1-7 
L 
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Extraction of acetic acid from water — 2 


Table 6 


Ethyl Acctate- Acetic Acid-W ater. 


Tie Line Data 


Ester Phase 

rc — 
% Acid % Ester | % Water % Acid 

30 3-0 90-1 6-0 3 

20-0 55-0 25-0 

9-6 80-4 10-0 0-6 

+0 89-1 6-0 5-2 

49 89-1 6-0 

0-8 80-3 9-9 10-0 

20-6 546 24°58 19-8 

3-90 90-1 6-0 

50 89-0 6-0 5-6 

10-1 80-4 5 10-0 

20°8 545 24-7 20-2 

10-0 10-0 80-0 | 9-8 


Water Phase Plait Point 
% Ester | %, Water | % Acid % Ester | % Water 

8-0 87-7 

23-0 

10-0 80-4 22:5 38-3 39-2 
8-0 86-8 

8-0 86-8 

10-0 80-0 

23-1 511 38-8 39-7 
8-0 87-7 

8-0 86-4 | 

10-0 80-0 

23-0 56-8 21-0 39-0 0-0 
10-0 80-2 


In the above table all percentages are by weight 


Prepiction oF Equitiprium Data 

The vapour liquid equilbria data of Brown 
and Ewa.p [21] for acetic acid-water and that 
of Pearce [22] for ethyl acetate-acetic acid, have 
been used to calculate ternary activity coefficients 
for acetic acid over a range of compositions on the 
solubility curve. The above binary activity 
coefficient plots were approximately corrected for 
temperature effects by assuming that T log y is 
constant. Since those binary systems are not 
satisfactorily correlated by the Margules and van 
Laar equations, linear interpolation has been 
used to calculate the ternary activity coefficients. 
The resulting activity composition plots are 
shown on Fig. 5. 


A comparison between predicted and experi- 
mental distribution ratios is given on Table 7. 

There is resonable agreement between the 
experimental and predicted results except for 
9-6 per cent by weight of acetic acid in the water 
phase. 

The above comparison between experimental 
and predicted results must, however, be treated 
with caution since the order of the activity 
coefficients for acetic acid in the two binary 
systems is quite low and thus susceptible to 
error. Furthermore there are association effects 
and some doubts as to the validity of using 
linear interpolation for predicting ternary activity 
coefficients. 


Table 7 


Acetic acid in water phase 


Weight %, Mole 
19-8 
9-6 3-4 
43 1-9 


Mole %, acid in ester phase 


Distribution Ratio = 


Mole % acid in water phase 


Experimental 


Predicted 
16 1°45 
2-9 2-1 
2-4 2-4 
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F. H. Garner and S. R. M. Ex.s : 


SOLUTROPY 


As the selectivity of acetic acid for ethyl! acetate 
reverses at a certain stage and then becomes 
selective towards water, Sounont and Waruap- 
PANDE [16] state this system is a solutrope at 
9-5 per cent by weight of acetic acid in the water 
layer. From Table 7 
molal basis the distribution coefficient is always 
greater than unity, and thus no molal solutrope 
exists, 

At 9-6 per cent by weight of acetic acid in the 


it can be seen that on a 


Extraction of acetic acid from water — 2 


water layer the mole per cent distribution ratio 
reaches a maximum value. Furthermore, at this 
point there is the greatest deviation between 
predicted and experimental distribution ratios. 
This discrepancy in the distribution ratios at the 
weight solutrope composition may imply there is 
a change on the degree of association in one of 
the phases. 

In the evaluation of the ternary activity 
coefficients the same molecular weight has been 
used for acetic acid in the two binary systems. 
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Book review 


Book review 


W.G. Bert: Physical Methods in Chemical Analysis, 
Academic Press, Inc. 640 pages. Price $13.50, 1951. 
As in Volume 1, each chapter in this book has been written 
by an expert in the field, and is concerned with a particular 
aspect of the use of physical methods in chemical analysis. 
The first chapter deals with polarography, and gives concise 
accounts of the theory, apparatus, and techniques used. 
Useful tables of depolarisation potentials are included, 
and over forty typical analyses are summarised. The 
chapter ends with a discussion of polarometric titrations. 
The basic features of polarography are adequately covered, 
and it is therefore unfortunate that only a few references 
to the available literature have been included for those 
who wish to follow the subject further. This is a general 
criticism of many chapters of the book, The next chapter 
covers the theory and technique of conductivity and 
conductometric titration. As a minor criticism one may 
note that the use of ion-exchange resins is not mentioned 
in the section on the preparation of conductivity water. 
The conductometric titration section includes phenols 
and alkaloids as well as inorganic acids and basis, and 
titrations which involve precipitations are also reviewed. 
Potentiometric analysis is discussed in an admirable 
manner in the following chapter, and adequate references 
to literature data and reference books are given. 


An excellent account of electrography and _ electro- 
spot testing includes a historical survey, full details, 
with good diagrams, of the equipment, materials, electro- 
lytes, reagents, etc. used and a good review of uses in 
identification of pure metal surfaces, alloys, and anions. 
Special applications, such as the structure of steel and non- 
ferrous metals, and a section on quantitative applications, 
complete this useful chapter. Adequate accounts of 
magnetic methods of analysis, and of the determination 


of the area of surfaces of solids, are followed by a chapter 
on surface tension measurements. This reviews com- 
petently the capillary rise, maximum bubble pressure, 
drop-weight and pendant drop, and ring methods of 
determination of surface tension, but would have been 
much improved by a detailed section on applications. 


Vacuum techniques of analysis depend first on the 
production and measurement of low pressures, and a good, 
but perhaps over-long, account of this is given. The 
vacuum fusion method for the analysis of gaseous elements 
in metals, and the determination of carbon by low-pressure 
combustion are next considered in detail. The section 
on the micro-analysis of gases at low pressures covers 
only a small part of the field, however, although literature 
references are given to other work. A discussion of 
molecular distillation and a full collection of references 
completes this somewhat unbalanced chapter. A good 
account of analysis by thermal conductivity comes next, 
and is followed by a useful chapter on the measurement 
of radioactivity for tracer applications. This surveys, 
with the aid of excellent diagrams, the photographic and 
electric methods of detection, the principles of counting, 
and analytical procedures. An equally useful chapter is 
that on statistical analysis, with particular reference to 
spectrography and spectrographic analysis; numerical 
examples are given. The final chapters surveys briefly 
chromatographic analysis, and is useful as an introduction 
to the subject. 

The book is well printed, is relatively free from major 
errors, and on the whole provides a valuable compendium 
of methods and techniques of physical methods of analysis 
which should be available in every library. 


R. N. HaszELDINE 
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Correction 


CORRECTION 


The condensation of vapour gas—vapour 
mixtures 


The constants to use in the equation on page 187 
of the June 1953 issue have been incorrectly 
evaluated. They should read : 


for chloroform and air K! — 6980 
for steam and air K' = 9190 


This introduces an error in the application of 
the Colburn approximate design method to the 
author's results. The new constants result in the 
calculated areas being approximately 90 per cent 
greater than the experimental area for both 
systems instead of approximately 40 per cent as 


reported. 
The conclusion that this method is conservative 
for the cases studied still applies. 
R. C. Carrns 
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